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ABSTRACT 

Swift captured for the first time a smoothly rising X-ray re-brightening of clear non- 
flaring origin after the steep decay in a long gamma-ray burst (GRB): GRB08f028. 
A rising phase is likely present in all GRBs but is usually hidden by the prompt tail 
emission and constitutes the first manifestation of what is later to give rise to the 
shallow decay phase. Contemporaneous optical observations reveal a rapid evolution 
of the injection frequency of a fast cooling synchrotron spectrum through the optical 
band, which disfavours the afterglow onset (start of the forward shock emission along 
our line of sight when the outflow is decelerated) as the origin of the observed re- 
brightening. We investigate alternative scenarios and find that the observations are 
consistent with the predictions for a narrow jet viewed off-axis. The high on-axis 
energy budget implied by this interpretation suggests different physical origins of the 
prompt and (late) afterglow emission. Strong spectral softening takes place from the 
prompt to the steep decay phase: we track the evolution of the spectral peak energy 
from the 7-rays to the X-rays and highlight the problems of the high latitude and 
adiabatic cooling interpretations. Notably, a softening of both the high and low spectral 
slopes with time is also observed. We discuss the low on-axis radiative efficiency of 
GRB 081 028 comparing its properties against a sample of Swift long GRBs with secure 
-E 7i i so measurements. 

Key words: gamma-ray: bursts - radiation mechanism: non-thermal -X-rays: indi- 
vidual (GRB081028). 



1 INTRODUCTION 

Gamma ray bursts (GRBs) are transient events able to out- 
shine the 7-ray sky for a few seconds to a few minutes. 
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The dis covery of their opt ical l|van Paradiis et al. .11 19971) and 
X-ray (|Costa et all 1 19971 ) long-lasting counterparts repre- 
sented a breakthrough for GRB science. Unfortunately, due 
to technological limitations, the X-ray observations were 
able to track the afterglow evolution starting hours after 
the trigger: only after the launch of the Swift satellite in 
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2004 (jGehrels et all |2004| ) was this gap between the end 
of the prompt emission and several hours after the onset 
of the explosion filled with X-ray obse rvations. A canonica l 
picture was then established (see e.g.. iNousek et al"1 [2006). 
with four different stages describing the overall structure 
of the X-ray afterglows: an initial steep decay, a shallow- 
decay phase, a normal decay and a jet-like decay stage. Er- 
ratic flares are found to be superimposed mainly to the first 
and second stage of emission. An interesting possibility is 
that the four light-curve phases instead be long to only two 
differ ent components of emission (see e.g., IWillingale et al.l 
120071 ): the first, connected to the activity of the central en- 
gine giving rise to the p rompt emission, comprises the flares 
l|Chincarini et all2007l and references therein) and the steep- 
decay phase; the second is instead related to the interaction 
of the outflow with the external medium and manifests it- 
self in the X-ray regime through the shallow, normal and 
jet-like decay. Observations able to further characterise the 
two components are therefore of particular interest. 

The smooth connection of the X-ray steep decay light- 
curve phase with the prompt 7-ray emission strongly sug- 
gests a common p hysical origin (jTagliaferri et al.l 1200a : 
lO'Brien et~aH l200fj): the hi g h latitude emission (HLE) 
model ( Fenimore et al.l 1 19961 : iKumar fc Panaitescul |2000| ) 
predicts that steep decay photons originate from the de- 
lay in the arrival time of prompt emission photons due to 
the longer path length from larger angles relative to our 
line of sight, giving rise to the a — j3 + 2 relation (where 
a is the light-curve decay index and (5 is the spectral en- 
ergy index). No spectral evolution is expected in the sim- 
plest formulation of the HLE effect in the case of a sim- 
ple power-law prompt spectrum. Observations say the op- 
posite: significant variations of the photon index have been 
found in the m ajority of GRBs during the steep decay 
phase (see e.g., IZhang. Liang fc Zhandl2007l ); more than 
this, the absorbed simple power-law (SPL) has proved to 
be a poor description of the spectral energy distribution of 
the steep decay phase for GRBs with the best statistic^]. 
A careful analysis of these events has shown t heir spectra 
to be best fit by an evolving Band function ijBand et al.l 
1 19931 ). establishing the link between steep decay and prompt 
emission photons also from the spectral point of view 
(see e.g., GRB 0606 14, iMangano et all 120071 : GRB 070616, 
IStarling et al.1 12008:): caused by the shift of the Band spec- 
trum, a temporal steep decay p hase and a spectral softening 
appear simultan eously (see e.g. IZhang. Liang fc WandlioCM 
lOin et all [2009). In particular, the peak energy of the vF„ 
spectrum is found to evolve to lower values, from the 7-ray 
to the soft X-ray energy range. Both the low (as observed 
for GRB 070616) and the high-energy portion of the spec- 
trum are likely to soften with time, but no observation is 
reported to confirm the high energy index behaviour during 
the prompt and steep decay phase. The observed spectral 
evolution with time is an invaluable footprint of the physi- 
cal mechanisms at work: observations able to constrain the 



behaviour of the spectral parameters with time are therefore 
of primary importance. 

By contrast, no spectral evolution is observed 
in t he X-ray during the shal low decay phase (see 
e.g. iLiang. Zhang fc Zhand I2007T ) experienced by most 
GRBs between ~ 10 2 s and 10 3 - 10 4 s. An unex- 
pected discovery of the Swift mission, the shallow de- 
cay is the first light-curve phase linked to the sec- 
ond emission component. A variety of theoretical ex- 
planations have been put forwar d. The proposed mod- 
els include: en ergy injection llPanaitescu et all 12006 



Rees fc Meszaros 1998; Granot fc Kumar 2006; Zhang et all 



2006 1 ) : reverse shock (see e.g.. lGenet. Daigne fc Mochkovitch 
2007); time dependent micro-ph y sical parameters (see 



e.g. iGranot. Konigl fc Piranl l2006l: lloka, et all 120061 ): off- 
axis emission (|Eichler fc Granotl l2006h : dust scattering 
l|Shao fc Daill2007l ). The predictions of all these models can 
only be compared to observations tracking the flat and de- 
cay phase of the second emission component, since its rise 
is usually missed in the X-ray regime, being hidden by the 
tail of the prompt emission. 

GRB 081028 is the first and unique event for which Swift 
was able to capture the rise of the second emission compo- 
nent the time properties of its rising phase can be con- 
strained for the first time while contemporaneous optical 
observations allow us to track the evolution of a break en- 
ergy of the spectrum through the optical band. GRB 081028 
is also one of the lucky cases showing a spectrally evolv- 
ing prompt emission where the evolution of the spectral pa- 
rameters can be studied from 7-rays to X-rays, from the 
trigger time to ~ 1000 s. A hard to soft spectral evolution 
is clearly taking place beginning with the prompt emission 
and extending to the steep decay ph ase, as already found 
for other Sw ift GRBs (GRB 0606 14, IMangano et~atl 120071 : 
GRB 070616, [starling et ail 120081 . are showcases in this re- 
spect). Notably , for GRB 081028 a softening of the slope of a 
Band function ((Band et al.ll 19931 ) above E p is also observed. 

The paper is organised as follows: Swift and ground- 
based observations are described in Sect. [2] data reduction 
and preliminary analysis are reported in Sect. [3] while in 
Sect.[4]the results of a detailed spectral and temporal multi- 
wavelength analysis are outlined and discussed in Sect. [5] 
Conclusions are drawn in Sect.|S] 

The phenomenology of the burst is presented in the 
observer frame unless otherwise stated. The convention 
F v (y, t) oc u~ l3 t~ a is followed, where f3 is the spectral 
energy index, related to the spectral photon index F by 
r = P + 1. All the quoted uncertainties are given at 68% 
confidence level (c.L): a warning is added if it is not the 
case. The convention Q x = Q/10 x has been adopted in 
cgs units unless otherwise stated. Standard cosmological 
quantities have been adopted: Ho = 70kms~ 1 Mpc~ , 
fi A = 0.7, Q M = 0.3. 



1 The limited 0.3 — lOkeV spectral co verage of the Swift X-Ray 
Telescope, XRT l|Burrows et al]|2005h . and the degeneracy be- 
tween the variables of the spectral fit can in principle lead to 
the identification of an SPL behaviour in intrinsically non-SPL 
spectra with poor statistics. 



2 There are a handful of long GRBs detected by Swift with a pos- 
sible X-ray rise of non -flaring origin. Amo ng them: GRB 070328 , 
iMarkwardt et al.l ||2007|): GRB 08022 9A, ICannizzo et al] feOOSft ; 
GRB 080307. iPage et al] J2009f) (see iPage et alj|2009l and refer- 
ences therein). However, in none of these cases has an X-ray steep 
decay been observed. A smooth rise in the X-rays has been ob- 
served in the short GRB 050724. 
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2 OBSERVATIONS 

GRB 81028 triggered the Swift Burst Alert Telescope 
fBAT: IBarthelmv et al.l l 2005h on 2008-10-28 at 00:25:00 UT 
( Gui dorzi et al .1120081 ). The spacecraft immedi ately slewed to 
the b urst allowing the X-ray Telescope fXRT: lBurrows et al.l 
l2005h to collect photons starting at T + 191 s after the trig- 
ger: a bright and fading X-ray a fterglow was discover ed. The 
UV/Optical Telescope fUVOT. iRoming et alj 12005) began 
observing at T + 210 s. In the first orbit of observations, 
no afterglow candidate was detected in any of the UVOT 
filters in either the individual or co-added exposures. A 
careful re-analysis of the acquired data revealed the pres- 
ence of a source with a White band magnitude of 20.9 at 
~ T+270 s (this paper) . A refined position was quickly avail- 
able thanks to the XRT-UVOT alignment procedure and the 
match of UVOT fi eld sources to the USNO-B1 catalogue (see 
iGoad et all 120071 for details): R.A.(J2000)=08 h 07 m 34.76 s , 
Dec.(J2000) =+02°18'29.8" wi th a 90% error radius of 
1.5 arcsec l|Evans et all 120081 ). Starting at ~ T + 9ks 
the X-ray light-curve shows a remarkable re-brightening 
jGuidorzi et al.ll2008bl ). see Fig. [2] this was later detected in 
ground-based near-infrared (NIR) and optical observations. 
Preliminary analysis r esults for this burst were reported in 
iGuidorzi et"aH (|2008d\ 

The Telesco pe a Action Rapid e pour les Objets Transi- 
toires (TAROT: Iklotz et all [20081 1 began observing 566.4 s 
after the trigger under poor weather conditions: no variable 
source was detected down to R ~ 17.4. 

The optical afterglow was discovered by the Gamma- 
Ray Burst Optical and Near-Infrared Detector (GROND; 
iGreiner et al .1120081 ). The observations started 20.9 ks after 
the trigger: the after glow was simultaneo usly detected in the 
g'r'i'z'JHK bands ^Clemens et al.ll2008h with the following 
preliminary magnitudes: g = 19.9±0.1; r' — 19.3±0.1; i' = 
19.2 ±0.1; z = 19.1±0.1; J = 19.0±0.15; H = 18.7±0.15; 
K = 19. 0± 0.15, with a net exposure of 264 and 240 s for the 
g'r'i'z' and the JHK bands respectively. Further G ROND 
observations were reported bv lClemens et all l|2008bl ) 113 ks 
after the trigger with 460 s of total exposures in g'r'i'z' and 
480 s in JHK. Preliminary magnitudes are reported below: 
g = 21.26 ± 0.05; r' = 20.49 ± 0.05; i' = 20.24 ± 0.05; 
z' = 19.99 ±0.05; J = 19.6 ±0.1. The source showed a clear 
fading with respect to the first epoch, confirming its nature 
as a GRB afterglow. 

The Nordic Optical Telescope (NOT) imaged the field 
of GRB 081028 ~ 6 hr after the trigger and independently 
confirmed the optica l afterglow with a magnitude R ~ 19.2 
l|Qlofsson et al.l 12008). Because of the very poor sky condi- 
tions only 519 frames out of 9000 could be used, with a total 
exposure of 51.9 s. The average time for the observations is 
estimated to be 05:53:00 UT. Image reduction was carried 
out by following standard procedures. 

An UV/optical re- brightening was discovered by the 
UVOT starting T + 10 ks, simultaneous to the X-ray re- 
brightening. T he afterglow was detected in the v, b and 
u-band filters (|Shadv et all 12008? ). The UVOT photometric 
data-set of GRB 81028 is reported in Tab. |Bl] We refer to 
iPoole et all (|2008l ) for a detailed description of the UVOT 
photometric system. 

The rising optical afterglow was independently con- 
firmed by the Crimean telescope for Asteroid Observations 



0.06 




-100 100 200 300 400 

Observed time from trigger [s] 



Figure 1. Top panel: BAT 15-150 keV mask weighted light-curve 
(binning time of 4.096 s). Solid blue line: 15 — 150 keV light-curve 
best fit using Norris et al. (2005) profiles. The typical 1 a error 
size is also shown. Bottom panel: best fit photon index T 7 as a 
function of time (errors are provided at the 90% c.L). 

(CrAO) and by the P eters Automated Infrared Imaging 
Telescope (PAIRITEL; iBloom et ai"1l2006l ). CrAO observa- 
tions were carried out starting at ~ T ± 1 ks and revealed 
a sharp rising optical afterglow peaking after T + 9.4 ks: 
R = 21.62 ± 0.07 at t = T + 1.8 ks; I = 21.32 ± 0.09 
at t = T + 3.6 ks; / = 21.43 ± 0.09 at t = T + 5.5 ks; 
I = 21.20 ± .08 at t = T + 7.5 ks; I = 20.66 ± 0.05 at 
t = T±9.4ks (jRumvantsev et all [20081 ). 

PAIRITEL observations were carried out 40 ks after 
the trigger: the afterglow was simultaneously detected in 
the J, H, and K s filters with a preliminary photometry 
J = 17.7 ± 0.1, H = 17.0 ± 0.1 and K s = 16.1 ± 0.1 
|Miller et al.ll2008l ). A total of 472 individual 7.8 s exposures 
were obtained under bad conditions (seeing > 3") for a total 
exposure time of ~3682 s. The data were redu ced and anal- 
ysed using the standard PAIRITEL pipeline (IBloom et alj 
2006). Photometry calibration was done against the 2MASS 
system. The resulting fluxes an d magnitude s are consis- 
tent with the values reported by I Miller et a l. (2008): how- 
ever, this work should be considered to supersede the pre- 
vious findings. The ground-based photometric data-set of 
GRB 081028 is reported in Tab. [Bl while the photometric 
optical observations of GRB 081028 are portrayed in Fig. [2] 

A spectrum of the GRB 081028 afterglow was taken 
with the Magellan Echellette Spectrograph (MagE) on the 
Magellan/Clay 6.5-m telescope at ~ T ± 27 ks for a total 
integration time of 1.8 ks. The identification of absorption 
features including SII, NV, SilV, CIV and Fell allowed the 
measurement of the redshift z — 3.038 t ogether with the dis - 
covery of several intervening absorbers |Berger et al.ll2008l ). 

According to lSchlegel et alj (|1998) the Galactic redden- 
ing along the line of sight of GRB 081028 is E(B-V) = 0.03. 



3 SWIFT DATA REDUCTION AND 
PRELIMINARY ANALYSIS 

The BAT data have been processed using standard Swift- 
BAT analysis tools within heasof t (v. 6. 6.1). The gr ound - 
refined coordinates provided by IBarthelmv et al.l (|2008l ) 
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have been adopted in the following analysis. Standard fil- 
tering and screening criteria have been applied. The mask- 
weighted background subtracted 15 — 150 keV is shown in 
Fig. [1] top panel. The mask-weighting procedure is also ap- 
plied to produce weighted, background subtracted counts 
spectra. Before fitting the spectra, we group the energy 
channels requiring a 3-er threshold on each group; the thresh- 
old has been lowered to 2-cr for spectra with poor statistics. 
The spectra are fit within Xspec v. 12. 5 with a simple power- 
law with pegged normalisation (pegpwrlw). The best fit 
photon indices resulting from this procedure are shown in 
Fig. [T] bottom panel. 

XRT data have also been processed with HEASOFT (v. 
6.6.1) and corresponding calibration files: standard filtering 
and screening criteria have been applied. The first orbit data 
were acquired entirely in WT mode reaching a maximum 
count rate ~ 140 counts s _1 . We appl y standard pile- u p cor - 
rections following the prescriptions of I Romano et all (2006) 
when necessary. Starting from ~ 10 ks Swift-~XRT switched 
to PC mode to follow the fading of the source: events are 
then extracted using different region shapes and sizes in or- 
der to maximize the signal-to-noise (SN) ratio. The back- 
ground is estimated from a source free portion of the sky. 
The resulting X-ray light-curve is shown in Fig. [2] the dis- 
played data binning assures a minimum SN equals to 4 (10) 
for PC (WT) data. In this way the strong variability of WT 
data can be fully appreciated without losing information 
on the late time behaviour. We perform automatic time re- 
solved spectral analysis, accumulating signal over time inter- 
vals defined to contain a minimum of ~ 2000 photons each. 
The spectral channels have been grouped to provide a min- 
imum of 20 counts per bin. The Galactic column density in 



(weighted average value from the lKalberla et aL I l2005l map). 
Spectral fitting is done within Xspec (v. 12.5) using a photo- 
electrically absorbed simple power law (SPL) model. The 
Galactic absorption component is frozen at the Galactic 
value together with the redshift, while we leave the intrin- 
sic column density free to vary during the first run of the 
program. A count-to-fiux conversion factor is worked out 
from the best fit model for each time interval for which 
we are able to extract a spectrum. This value is consid- 
ered reliable if the respective xVdof (chi-square over de- 
grees of freedom) implies a P-value (probability of obtain- 
ing a result at least as extreme as the one that is actually 
observed) higher than 5%. The discrete set of reliable count- 
to-flux conversion factors is then used to produce a contin- 
uous count-to-flux conversion factor through interpolation. 
This procedure produces flux and luminosity light-curves 
where the possible spectral evolution of the source is prop- 
erly taken into account (Fig. [2}. In the case of GRB 081028 
this is particularly important: the simple power law photon 
index evolves from F~1.2tor~3 during the steep decay 
phase (Fig. |S} , inducing a variation of a multiplicative fac- 
tor ~ 1.7 in the count-to-flux conversion factor. As a second 
run, we remove one degree of freedom from the spectral fit- 
ting procedure, noting the absence of spectral evolution dur- 
ing the X-ray re-brightening in the X-ray regime (see Sect. 
14.40 . This gives the possibility to obtain a reliable estimate 
of the intrinsic neutral Hydrogen column density Nn,z of 
GRB 081028: the PC spectrum accumulated over the time 
interval 10 — 652 ks can be adequately fit by an absorbed 
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Table X. Best fit parameters and related quantities resulting from 
the modelling of the prompt 15-150 keV emission with two Norris 
et al. (2005) profiles. From top to bottom: peak time, start time, 
1/e rise time, 1/e decay time, 1/e pulse width, pulse asymmetry, 
peak count-rate and statistical information. The x 2 value mainly 
reflects the partial failure of the fitting function to adequately 
model the peaks of the pulses (see Norris et al. 2005 for details). 



SPL model with best fit photon index V = 2.09 ± 0.07 and 
Nn,* = (0.52 ± 0.25) x 10 22 cm -2 (90% c.l. uncertainties 
are provided). The flux-luminosity calibration procedure is 
then re-run freezing the intrinsic absorption component to 
this value. 

T he UVOT photom etry was performed using standard 
tools (|Poole et alJl200Sl ) and is detailed in Tab. [BT] 



4 ANALYSIS AND RESULTS 

4.1 Temporal analysis of BAT (15-150 keV) data 

The mask-weighted light-curve consists of two main pulses 
peaking at T + 70 s and ~ T + 200 s followed by a long 
lasting tail out to ~ T + 400 s. In the time interval T — 100 s 
T±400s, th e light -curve can be fit by a combination of two 
iNorris et al.l (|2005l ) profiles (Fig. [T] top panel) , each profile 
consisting of the inverse of the product of two exponentials, 
one increasing and one decreasing with time. The best fit 
parameters and related quantities are reported in Tab le [1] 
the parameters are defined following INorris et all (|2005l ): we 
account for the entire covariance matrix during the error 
propagation procedure. The GRB prompt signal has a T90 
duration of 261.0 ± 28.7 s and a T 50 = 128.2 ± 7.7 s. 

The temporal variability of this burst ha s been charac- 
terised in two different ways. First, following iRizzuto et al.l 
(2007) we compute a variability measure Var(15 — 
150 keV) = (5.0 ± 0.14) x 10 -2 . Second, we ado pt the 
power spectru m analysis in the time domain (|lT l200ll : 
iLi et alj 120021 ) : unlike the Fourier spectrum, this is suit- 
able to study the rms variations at different time-scales. See 
iMargutti et alj (|2008h and Margutti et al. in prep, for details 
about the application of this technique to the GRB prompt 
emission. In particular, we define the fractional power den- 
sity (fpd) as the ratio between the temporal power of the 
source signal and the mean count rate squared. This quan- 
tity is demonstrated to show a peak at the characteristic 
time scales of variability of the signal. We assess the signifi- 
cance of each fpd peak via Montecarlo simulations. The fpd 
of GRB 081028 shows a clear p eak around 70 s (t ime scale 
related to the width of the two INorris et al.ll2005l profiles). 
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Figure 2. Complete data set for GRB 081028 starting 200 s after the trigger including X-ray (XRT, flux density estimated at 1 keV), 
UV/visible/NIR (UVOT, GROND, PAIRITEL, CrAO, NOT) observations. The arrows indicate 3-cr upper limits of UVOT observations. 
The shaded regions indicate the time intervals of extraction of the SEDs. 



Below 70 s the fpd shows a first peak at At ~ 2 s and then 
a second peak at At ~ 6 s, both at 1-a c.l. The signal shows 
power in excess of the noise at 2-cr c.l. significance for time 
scales At > 32 s. 



4.2 Spectral analysis of BAT (15-150 keV) data 

We extract several spectra in different time intervals and 
then fit the data using different models to better constrain 
the spectral evolution of GRB 081028 in the 15-150 keV en- 
ergy band. The first spectrum is extracted during the Tgo 
duration of the burst; a second spectrum is accumulated dur- 
ing the entire duration of the 15-150 keV emission; finally, 
the signal between 10 s and 290 s from trigger has been split 
into two parts, taking 150 s as dividing time, to characterise 
the spectral properties of the two prompt emission pulses. 
The resulting spectra are then fit using a simple power-law 
and a cut-off power-law models within Xspec. The results 
are reported in Table [2] The measured simple power law 
photon index around 2 suggests that B AT observed a por - 
tion of an intrinsically Band spectrum dBand et alj|l993l ). 
Consistent with this scenario, the cut-off power law model 
always provides a better fit which is able to constrain the 
peak energy value (E p , peak of the vF v spectrum) within 
the BAT energy range. 

The best fit parameters of the cut-off power-law model 
applied to the total spectrum of Table [2] imply i?i BO , 7 (l — 



10 4 keV) = (1.1 ± 0.1) x 10 53 erg. The respective rest frame 
peak energy is E p> i — (1 + z)E p = 222l|j!j keV, placing 
GRB 081028 within the 2-a region of the Amati relation 
The burst is characterised by an isotropic 
10 2 -10 3 keV (rest frame) i iso = (2.85T0.25) x 10 51 ergs" 1 . 
This information together with the variability measure 
Var(15- 150 keV) = (5.0 ±0.14) x 10" 2 makes GRB 081028 
per fectly consistent with the luminosity variability relation 
(see l Reichart et al.ll200ll ; iGuidorzi et al.ll2005l ; iRizzuto et al.l 
l2007| y 



4.3 Temporal analysis of XRT (0.3-10 keV) data 

The XRT (0.3-10 keV) light-curve consists of two parts: a 
steep decay phase with flares and variability superimposed 
(100 s < t < 7000 s), followed by a remarkable re-brightening 
with smoothly rising and decaying emission between 7 ks 
and 1000 ks. The two light-curve phases are studied sepa- 
rately. 

GRB 081028 is one of the rare cases in which the XRT 
caught the prompt emission. The light-curve shows a flat 
phase up to t ~ 300 s followed by a steep decay. Starting 
from ~ 690 s the light-curve is dominated by a flare which 
peaks at 800 s but whose decaying phase is temporally coin- 
cident with the orbital data gap. The steep decay behaviour 
before the flare is inconsistent with the back-extrapolation 
of the post orbital data gap power-law decay, as shown in 
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Interval 


Model 


^start 

00 


^stop 

00 


V,a 


E P 
(keV) 


Fluencc 
(erg cm -2 ) 


X 2 /dof 


P-value 


T'M) 


PI 

Cutpl 


52.9 
52.9 


317.2 
317.2 


1.82 ± 0.09 
1.3 ±0.4 


65+« 


(3.3 ± 0.20) X 10~ 6 
(3.15 ± 0.20) x IO" 6 


31.8/31 
25.8/30 


43% 
69% 


Total 


PI 
Cutpl 


0.0 
0.0 


400.0 
400.0 


1.89 ±0.09 
1.3 ±0.4 


55±f 


(3.7 ±0.20) x 10" 6 
(3.45 ± 0.19) x 10~ 6 


37.4/32 
30.1/31 


23% 
51% 


Pulse 1 


PI 
Cutpl 


10.0 
10.0 


150.0 
150.0 


1.91 ±0.13 
1.1 ±0.6 


49t^ 8 


(1.60 ± 0.12) x 10" 6 
(1.47 ±0.11) x IO" 6 


18.0/24 
12.0/23 


80% 
97% 


Pulse 2 


PI 
Cutpl 


150.0 
150.0 


290.0 
290.0 


1.77±0.11 
1.22 ±0.45 


69ifl 


(1.79 ±0.11) x 10~ 6 
(1.47 ±0.11) x 10~ 6 


33.8/29 
29.3/28 


25% 
40% 



Table 2. Best fit parameters derived from the spectral modelling of 15-150 keV data using a power law with pegged normalisation (PI, 
PEGPWHLW within Xspec) and a cut-off power-law model with the peak energy of the vF v spectrum as free parameter (Cutpl). From left 
to right: name of the interval of the extraction of the spectrum we refer to throughout the paper; spectral model used; start and stop 
times of extraction of the spectrum; best fit photon index Y for a PI model or cutoff power-law index for a Cutpl model; best fit peak 
energy of the vF v spectrum; fluence; statistical information about the fit. 





C 


m 


a 


b 


di 


<bri 

(ks) 


"3 


e d 2 


^br 2 

(ks) 


X 2 /dof 


-ml8.3±5.5 
1 g29.9±6.5 


-5.2 ± 1.5 
-7.60 ± 1.8 


1.2 ± 1.1 
0.31 ±0.02 


-4.5 ±3.3 
-1.8 ±0.3 


2.1 ±0.1 
1.3 ±0.1 


2.4 ± 2.0 
0.1 


15.5 ±6.3 
19.5 ±0.7 


0.06 


2.3 ±0.1 0.05 


62 


164.8/145 
147.1/143 



Table 3. Best-fit parameters of the XRT light-curve modelling starting from 3 ks after the trigger. The first (second) line refers to Eq. 
Q](Eq.[2}. 
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Figure 3. 0.3-10 keV X-ray afterglow split into different compo- 
nents. Green dot-dashed line: steep decay; purple long dashed 
line: pre-rebrightening component; light grey region: first re- 
brightening component; dark grey region: second re-brightening 
component; red solid line: best fit model. See Sect. |4~3.1I for de- 
tails. 



Fig. [3] The strong spectral evolution detected by the XRT 
(Sect. 14. 4p requires a time resolved flux calibration of the 
light-curve before the light-curve fitting procedure. In the 
time interval 320 s < t < 685 s the 0.3-10 keV light-curve 
best fit is given by a simple power- law with a — 3.6 ±0.1 
(x 2 /dof = 768.3/736). Figure g] shows the different tem- 
poral behaviour of the detected signal when split into dif- 
ferent energy bands: harder photons decay faster. The 0.3- 
1 keV light-curve decays following a power-law with index 
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Figure 4. Upper panel: steep decay portion of GRB 081028 X- 
ray afterglow. The XRT signal has been split into 3 energy bands 
so that the different temporal behaviour can be fully appreciated. 
Lower panel: (0.3 — lkeV)/(l — lOkeV) hardness ratio evolution 
with time. The signal clearly softens with time. In both panels, 
the vertical black dashed lines mark the orbital data gap. 

q w —2.5; the decay steepens to a w —3.5 and a w —3.8 for 
the 1-2 keV and 2-10 keV signal, respectively. 

During the re-brightening there is no evidence for spec- 
tral evolution in the XRT energy band (see Sect. 14. A\ . For 
this reason we model the count-rate light-curve instead of 
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Log[Obs Time /s] 

Figure 5. XRT 0.3-10 koV count-rate light-curve of GRB 081028 
starting from 3 ks with best fit model superimposed (Eq. [2] Table 
rSl . Inset: residuals with respect to the best-fit model. 



the flux calibrated one: this gives the possibility to obtain 
a fully representative set of best fit parameter^ determined 
with the highest level of precision. The count-to-flux cali- 
bration introduces additional uncertainty inherited by the 
spectral fitting procedure. Starting from 3 ks (the inclusion 
of the last part of the steep decay is necessary to model the 
rising part of the re-brightening) , the count-rate light-curve 
can be modelled by a p ower-law plus Beuermann function 
|Beuermann et aTlll999l ) where the smoothing parameter di 
is left free to vary: 



n 2 ■ t + m 



+ 



(1) 



t 

. tbri / \ tbr 

The best fit parameters are reported in Table [3] The draw- 
back of this model is that the best-fit slopes are asymptotic 
values and do not represent the actual power-law slopes. 
While due to the smooth transition between the rising and 
decaying phases, this makes the comparison between obser- 
vations and model predictions difficult. Freezing di at 0.1 
to have a sharp transition results in an unacceptable fit (P- 
value~ 10 -4 ) and suggests a light-curve steepening around 
50 ks. The possibility of a break is investigated as follows: 
we select data points starting from 20 ks and fit the data us- 
ing a SPL or a broken power-law (BPL) model. Given that 
the SPL and BPL models are nested models and the possi- 
ble values of the second model do not lie on the boundary 
of de finition of the first model parameters l|Protassov et al.l 
2002), we can apply an F-test: with a probability of chance 
improvement ~ 1%, we find moderate statistical evidence 
for a break in the light-curve at 62 ks. The final fitting func- 
tion is given by Eq. 



n 2 • t c + m 
/ • n 3 



I; 

d- 2 



+ 



-f-1 -d2 

d 2 



t < 40 ks 



t > 40 ks 



(2) 



where / is function of the other fitting variables and assures 
the continuity of the fitting function at 40 ks. The light-curve 
of GRB 081028 fits in this case with xV dof = 147.1/143 



3 This is in general not true in cases of strong spectral evolution 
as shown in the first part of this section. 



and an P-value=39%: the best fit parameters are reported 
in Table [3] while a plot of the result is provided in Fig. [5] 
The fit of the flux-calibrated light-curve gives completely 
consistent results. The model predicts -Fx, P = (1.53±0.08) x 
10 -11 ergcm~ 2 s _1 , where Fx, P is the flux at the peak of the 
re-brightening. 



4-3.1 Count-rate drop around 250 ks 

The drop of the count-rate around 250 ks is worth attention: 
the statistical significance of this drop is discussed below. We 
select data with t > 60 ks. These data can be fit by a sim- 
ple power-law with index a = 1.9 ± 0.2 (x 2 /dof = 11.0/12, 
P-value=53%). According to this model the drop is not sta- 
tistically significant (single trial significance of ~ 2.6 a). 
However, this model under-predicts the observed rate for 
t < 60 ks: an abrupt drop of the count-rate during the or- 
bital gap at 80 ks would be required in this case. Alterna- 
tively there is not any kind of switch-off of the source during 
the orbital gap and the flux around 80 ks joins smoothly to 
the flux component at t < 60 ks, as portrayed in Fig. [5] A 
careful inspection of the figure reveals the presence of a non 
random distribution of the residuals of the last 14 points, 
with the points before 250 ks being systematically low and 
those after 250 ks being systematically high. While this fit 
is completely acceptable from the \ 2 point of view, a runs 
test shows that the chance probability of this configuration 
of residuals is less than 0.1%. This would call for the intro- 
duction of a new component to model the partial switch-off 
and re-brightening of the source around 250 ks. A possi- 
ble description of the light-curve behaviour for t > 20 ks 
(peak time of the main re-brightening) is represented by a 
Beuermann plus Beuermann function with smoothing pa- 
rameters frozen to give sharp transitions; the first decay- 
ing power-law index is frozen to 6 = 1.3 ± 1.3 while the 
break time of the first Beuermann component is frozen to 
tbr 2 = 62 ks as reported in Table [3] The light-curve decays 
with a 2 = 3.1 ± 0.2 (a 3 = 1.5 ± 0.7) for 60 ks < t < 250 ks 
(t > 316 ks), see Fig. This additional component would 
account for ~ 10% of the total re-brightening 0.3 — 10 keV 
energy which is ~ 1.1 X 10 52 erg. 

The temporal properties of the secon d re-brightening 
seem to point to refreshed shocks (s ee e.g. iKumar fc Piranl 
l2000bl : iGranot. Nakar fc Piranl 120031 ): the decaying power- 
laws before and after the drop are roughly consistent with 
each other but shifted upwards in the count-rate axis. Since 
at this epoch the observed X-ray frequencies are above 
both the cooling and the injection frequencies, in the stan- 
dard afterglow scenario the X-ray flux is oc E^ 2 ^ i in- 
dependent of the external medium density profile (see e.g. 
iPanaitescu fc Kumar|[2000l . their appendices B and C): the 
observed jump in flux would therefore require an increase of 
the energy in the forward shock by a factor of ~ 3. Given 
the marginal statistical evidence, the properties of the sec- 
ond re-brightening will not be discussed further. 



4.4 Spectral analysis of XRT (0.3-10 keV) data 

The very good statistics characterising the X-ray afterglow 
of GRB 081028 gives us the possibility to perform a tem- 
porally resolved spectral analysis. Figure [6] shows the dra- 
matic evolution of the photo-electrically absorbed simple 



8 Margutti et al. 





















i 






O 

o 

o 


? % - 









10 3 10 4 10 5 10 6 

Obs Time (s) 



Figure 6. 0.3-10 keV light-curve (grey points, arbitrary units) 
with best fit 0.3-10 keV photon index superimposed (black 
points). Each point comes from the fit of a spectrum consisting 
of ~ 2000 photons: the model tbabs*ztbabs*pow within Xspec 
with the intrinsic column density JVh z frozen to 0.52 X 10 22 cm -2 
has been used. An exception is represented by the first data point 
after the orbital gap: see Sect. 14.41 for details. The vertical red 
dashed lines mark the time interval of the first orbital gap. An 
abrupt change of the spectral properties of the source temporally 
coincident with the onset of the re-brightening is apparent. 



power-law photon index with time during the first 1000 s 
of observation, with F evolving from 1.2 to 2.7. The intrin- 
sic neutral Hydrogen column density JVh.z has been frozen 
to 0.52 x 10 22 cm -2 for the reasons explained below. If left 
free to vary, this parameter shows an unphysical rising and 
decaying behaviour between 200 s and 600 s. 

The temporal behaviour of the light-curve in the time 
interval 4—7.5 ks (after the orbital gap, see Fig. [3} physically 
connects these data points with the steep decay phase. We 
test this link from the spectroscopic point of view. The 0.3- 
10 keV spectrum extracted in this time interval contains 
133 photons. Spectral channels have been grouped so as to 
have 5 co unts per bin and then weighted using the Churazov 
method (jChurazov et al.|[l996T ) within Xspec. A fit with a 
photo-electrically absorbed power-law (tbabs*ztbabs*pow 
model) gives V = 2.63 ± 0.25, (90% c.l., xV d °f = 25.6/23, 
P-value=32%), confirming that this is the tail of the steep 
decay detected before the orbital gap as shown by Fig. [II] 

The light-curve re-brightening around 7ks translates 
into an abrupt change of the 0.3-10 keV spectral prop- 
erties (Fig. 0, with T shifting from 2.7 to 2. The pos- 
sibility of a spectral evolution in the X-ray band during 
the re-brightening is investigated as follows: we extracted 
three spectra in the time intervals 7-19.5 ks (sped, rising 
phase); 19.5-62 ks (spec2, pre-break decaying phase); 62 ks- 
end of observations (spec3, post-break decaying phase). A 
joint fit of these spectra with an absorbed simple power- 
law model (tbabs*ztbabs*pow model) where the intrin- 
sic Hydrogen column density is frozen to 0.52 x 10 22 cm~ 2 
(see Sect. [3} and the photon index is tied to the same 
value, gives T = 2.04 ± 0.06 with xVdof = 118.0/167. 
Thawing the photon indices we obtain: Ti = 2.13t°'\t'> 
T 2 = 2.03±Ht; r 3 = 2.00± ;il (xVdof = 115.8/165). Un- 



certainties are quoted at 90% c.l.. The comparison of the 
two results implies a chance probability of improvement of 
22%: we conclude that there is no evidence for spectral evo- 
lution during the re-brightening in the 0.3-10 keV energy 
range. The same conclusion is reached from the study of the 
(1 - 10keV)/(0.3 - 1 keV) hardness ratio. 

4.5 Spectral energy distribution during the 
re-brightening: evolution of the break 
frequency 

The re-brightening properties can be constrained through 
the study of the temporal evolution of the spectral energy 
distribution (SED) from the optical to the X-ray. We extract 
4 SEDs, from the time intervals indicated by the shaded 
bands in Fig. [2] 

(i) SED 1 at t ~ 10 ks corresponds to the rising portion 
of the X-ray re-brightening and includes XRT and UVOT 
observations; 

(ii) SED 2 is extracted at t ~ 20 ks, peak of the X-ray 
re-brightening. It includes XRT, UVOT, GROND and NOT 
observations; 

(iii) SED 3 at f ~ 41 ks describes the afterglow spec- 
tral energy distribution during the decaying phase of the 
re-brightening, before the detected light-curve break. It in- 
cludes X-ray data from ~ 30 ks to ~ 62 ks, UVOT and 
PAIRITEL observations; 

(iv) SED 4 corresponds to the post-break decaying por- 
tion of the re-brightening, at t ~ 112 ks and includes XRT 
and GROND observations. 

When necessary, optical data have been interpolated to 
the time of extraction of the SED. Uncertainties have been 
propagated accordingly. 

At a redshift of 3.038, we expect some contamination in 
the spectrum from abso rption system s located between the 
Earth and GRB 081028 (|Madaulll995l ). This means that the 
g filter of GROND and all UVOT filters but the v band are 
marginally or strongly affected by Lyman absorption: these 
filters are consequently excluded from the following analysis. 

The Galactic and intrinsic absorption at wavelengths 
shorter than the Lyman e dge are modelled using th e photo- 
electric cross-sections of M orrison fc McCammonl (l983h . 
We ado pt t h e ana lytical description of the Galactic extinc- 
tion by [Pel <|l992h . while the host galaxy absorption is as- 
sume d to be mo delled by a Small Magellanic Cloud-like law 
(from |pei|[l992l ). 

An absorbed SPL model from the optical to the X-ray 
range is not able to account for SED 1, SED 2 and SED 3 
(Fig. 0, while it gave the best fit model for SED 4. For the 
first three SEDs a satisfactory fit is given by a broken power- 
law with X-ray spectral index /3 X ~ 1; optical spectral index 
Po ~ 0.5 and ATh,z consistent with the value reported in Sect. 
EH (0.52 x 10 22 cm" 2 ). The best fit break frequency is found 
to evolve with time to lower values following a power-law 
evolution with index a ~ 2. This evolution is faster than 
expected for t he cooling frequency of a synchrotro n spec- 
trum (see e.g. ISari et"aLlll998l : iGranot fc Sarill2002h : in the 
following, we identify the break frequency with the injec- 
tion frequency. We freeze the Galactic c ontribution to give 
E(B - V) = 0.03 (jSchleeel et al.lll998l ). while leaving the 
intrinsic component free to vary. 
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Figure 7. Observer-frame SED1, SED2, SED3 and SED4 from optical to X-ray extracted at t ~ 10 ks, t ~ 20 ks, t ~ 41 ks and t ~ 112 ks, 
respectively. Red solid line: photo-electrically absorbed model corresponding to Eq. \3\ This proved to be the best fit model for SED1, 
SED2 and SED3. Blue dashed line: photo-electrically absorbed simple power law. This is the best fit model for SED4. For all SEDs an 
SMC extinction curve at the redshift of the source is assumed. The best fit parameters are reported in Table [4] 



Th e broken power- l aw m odel has been then refined as 
follows. iGranot fc SarH (2002) showed that under the as- 
sumption of synchrotron emission from a relativistic blast 
wave that accelerates the electrons to a power law distri- 
bution of energies iV(7 e ) oc 7,T P , it is possible to derive 
a physically motivated shape of spectral breaks. Interpret- 
ing the break frequency as the injection frequency in the 
fast cooling regime, the broken power-law model reads (see 
IGranot fc SariH2002l , their Eq. 1): 

F »= F iU + U) \ (3) 

where Vb and F n are the break frequency and the normal- 
isation, respectively; /3i = —0.5 and @2 = —p/2 are the 
asymptotic spectral indices below and above the break 
under the conditions above; s = s(p) is the smoothing 
parameter: in particular, for an inters tellar (wind) medium 
s = 3.34 - 0.82p (s = 3.68 - 0.89p) (|Granot fc Sarill2002j . 
their Table 2). The free parameters of the final model are 
the following: normalization of the spectrum F n , break 
frequency v^, power- law index of the electron distribution 
p, intrinsic neutral Hydrogen column density Nh,z, and host 
reddening. The ISM or wind environments give perfectly 
consistent results. We choose to quote only ISM results for 
the sake of brevity. For SED 4 we use an absorbed simple 
power-law with spectral index — p/2. The four SEDs are 



SED 




Parameter 


Value 


1,2,3,4 




P 


1.97 ±0.03 


1 


Lof 


!l0 K/10 15 Hz) 


2.0 ±0.1 


2 


Lo{ 


!l0 K/10 15 Hz) 


1.4 ±0.1 


3 


Lo{ 


; 10 K/10 15 Hz) 


0.4 ±0.1 


X 2 /dof 
P-value 




134.7/138 
56% 





Table 4. Best fit parameters for the simultaneous fit of SED1, 
SED2, SED3 and SED4. For SED1, SED2 and SED3 the emission 
model is expressed by Eq. \3\ while for SED4 we used a simple 
power-law with spectral index p/2. The spectral normalisations 
and break frequencies have been left free to take different values 
in different spectra. The intrinsic neutral Hydrogen column value 
is found to be consistent with the value inferred from the X-ray 
spectra. 



first fit separately; as a second step we perform a joint fit 
where only the spectral normalisation and break frequency 
are free to take different values in different spectra. We find 
fully consistent results with improved uncertainties thanks 
to the tighter constraints imposed by the joint fit. The best 
fit results are reported in Table [4] and portrayed in Fig. [7] 
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Figure 8. Spectral break frequency (See Eq. |3)l evolution with 
time as found from a simultaneous fit of SED1, SED2, SED3 and 
SED4 with best fit models superimposed. Red solid line (blue 
dashed line): simple power-law with zero time to=0 ks (2 ks) and 
best power-law index a = 2.6 (2.3). The satisfactory fit of SED4 
with a simple power-law provides the upper limit shown. 

The spectral break frequency v\> evolves with time to 
lower values, as shown in Fig. [8] The consistency of SED4 
optical and X-ray data with a simple power-law model with 
index —p/2, suggests that the break frequency has crossed 
the optical band by the time of extraction of SED4. This 
translates into Log 10 (^ b /10 15 Hz) < -0.33 for t > 112 ks. 
The decrease of the break frequency with time can be mod- 
elled by a simple power-law function: this leads to an ac- 
ceptable fit (x 2 /dof = 1.4/1, P-value=24%) with best fit 
index a — 2.6 ± 0.2. Using to = 2 ks as zero time of the 
power-law model we obtain: a = 2.3 ± 0.2 (x 2 / dof = 2.1/1, 
P-value=15%). 

The fit implies a limited rest frame optical extinction 
which turns out to be E(B — V) z ~ 0.03. A 3 — a upper limit 
can be derived from the joint fit of the four SEDs, leaving all 
the parameters but the one related to the optical extinction 
free to vary. The upper limit is computed as the value which 
increases the \ 2 by a Ax 2 corresponding to a 3 a c.L. This 
procedure leads to: Av, z < 0.22. 

4.6 Peak energy evolution with time 

The consistency of the prompt BAT spectrum with a cut-off 
power-law (Sect. [5} and the spectral variability detected in 
the XRT energy range (Sect. 14.41 Fig. [SJ) suggests that the 
peak of the vF v spectrum is moving through the BAT+XRT 
bandpass. To follow the spectral evolution, we time slice the 
BAT and XRT data into 14 bins covering the 10-851 s time 
interval. The spectra are then fit within Xspec using a Band 
function (ngrbep) or a cut-off power-law (cutplep) with 
E p as free parameter; alternatively a simple power law is 
used. Each model is absorbed by a Galactic (hydrogen col- 
umn density frozen to 3.96 x 10 20 cm -2 ) and intrinsic compo- 
nent (A H ,z frozen to 0.52 x 10 22 cm" 2 , see Sect.gj}. When 
possible we take advantage of the simultaneous BAT and 
XRT observations, performing a joint BAT-XRT spectral 
fit. The normalisation for each instrument is always tied to 
the same value. The best fit parameters are reported in Ta- 
ble [S] the simple power law model gives a poor description 
of the spectra up to ~ 400 s, as the curvature of the spectra 
requires a cut-off power-law or a Band function. In partic- 
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Figure 9. Time resolved combined analysis of XRT and BAT 
data. Upper panel: BAT 15-150 keV and XRT 0.3-10 keV flux 
light-curves. No extrapolation of the BAT data into the XRT en- 
ergy range has been done. The vertical dashed lines mark the 
intervals of extraction of the spectra: these are numbered accord- 
ing to Table \E\ first column. Central panel: best fit photon in- 
dices evolution with time. Lower panel: best fit E p parameter as 
a function of time. The decay has been fit with a simple power- 
law model starting from 200 s from trigger: E p (t) oc (t — io) a • 
Starting from 405 s E p is likely to be outside the XRT energy 
range: E p < 0.3 keV (solid green line). 



ular, this is the case when the high energy slope enters the 
XRT bandpass. The E p parameter is well constrained and 
evolves to lower energies with time; at the same time both 
the high and low energy photon indices are observed to grad- 
ually vary, softening with time (Fig. [9} . The E p decay with 
time can be modelled by a simple power-law starting ~ 200 
s after trigger: E v oc (t — to) a . The best fit parameters are 
reported in Table [6] 

The uncertainty of the inter-calibration of the BAT and 
XRT has been investigated as possible source of the detected 
spectral evolution as follows. For each time slice, we multiply 
the fit model by a constant factor which is frozen to 1 for the 
BAT data. For XRT, this factor is left free to vary between 
0.9 and 1.1, conservatively allowing the XRT calibration to 
agree within 10% with the BAT calibration. The best fit 
parameters found in this way are completely consistent with 
the ones listed in Table [5] The inter-calibration is therefore 
unlikely to be the main source of the observed evolution. 
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Model 




08 (T) 


Bp 
(keV) 


X 2 /dof 


P-value 


3 


203 


222 


BAT+XRT 


Cutpl 
PI 


1.19 + 0.05 


1.37 + 0.02 


i on n 
61.01-2°;° 


100.4/112 
147.3/114 


77% 
2% 


4 


222 


247 


BAT+XRT 


Cutpl 
PI 


1.28 + 0.06 


1.44 ± 0.03 


41.5±^ 


80.5/101 
108.2/102 


93% 
31%+ 


5 


247 


271 


BAT+XRT 


Cutpl 
PI 


1.38 + 0.17 


1.54 + 0.04 


16.1111 


77.4/88 
96.5/89 


78% 
3%+ 


6 


271 


300 


BAT+XRT 


Cutpl 
PI 


1.57 + 0.07 


1.76 ± 0.03 


12.51" 


129.8/91 
158.6/92 


1% 
0.001% 


7 


300 


323 


XRT 


Cutpl 
PI 


1.20 + 0.16 


1.49 ± 0.05 


5 2+ 3 - 7 


77.1/81 
87.6/82 


60% 
32% 


8 


323 


343 


XRT 


Cutpl 
PI 


0.82 + 0.18 


1.61 ± 0.05 


2.9±g;i 


78.7/83 
149.8/84 


61% 
0.001% 


9 


343 


371 


XRT 


Cutpl 
PI 


1.38 + 0.17 


1.91 ± 0.05 


2.0+0| 

— U.o 


94.3/84 
131.2/82 


15% 
0.1% 


10 


371 


405 


XRT 


Band 
Cutpl 
PI 


~ 1.10 
1.81 ± 0.016 


o n+0.1 

z -°-0.2 

2.07 + 0.06 


< l.i 

-i ,-y+O 3 


82.4/77 
102.4/78 
109.7/79 


31% 
3% 
1% 


11 


405 


456 


XRT 


PI 




2.32 ± 0.06 




100.1/78 


5% 


12 


456 


530 


XRT 


PI 




2.34 + 0.06 




103.3/79 


3% 


13 


530 


664 


XRT 


PI 




2.61 ± 0.07 




98.1/76 


5% 


14 


664 


838 


XRT 


PI 




2.71 + 0.06 




89.3/73 


7% 


15 


838 


851 


XRT 


PI 




2.68 + 0.18 




15.7/10 


1% 



Table 5. Best fit parameters derived from the spectral modelling of XRT and BAT data using photo-electrically absorbed models 
(tbabs*ztbabs within Xspec). The BAT and XRT normalisations are always tied to the same value. Three different models have been 
used: a simple power-law (PI); a cut-off power law and a Band function both with the peak energy of the vF v spectrum as free parameter. 
From left to right: name of the interval of the extraction of the spectrum we refer to throughout the paper (intervals 1 and 2 correspond 
to Pulse 1 and Pulse 2 of Table|2j; start and stop time of extraction of each spectrum; energy range of the fit: "XRT+BAT" stands for 
a joint BAT-XRT data fitting; model used; best fit low and high energy photon indices for a Band or Cutpl power-law or best fit photon 
index T for a PI model; statistical information about the fit. The + symbol indicates an apparent trend in the residuals of the fit. 



(Chincarini et al., in prep.): only in GRB051016B is a ris- 
ing feature detected at the end of the steep decay Q. In this 
case, the sparseness of the data prevents us from drawing 
firm conclusions, so that a flare origin of the re-brightening 
cannot be excluded. 

The very good statistics of GRB 081028 allows us to 
track the detailed spectral evolution from 7-rays to X-rays, 
from the prompt to the steep decay phase: this analysis fully 
qualifies the steep decay as the tail of the prompt emis- 
sion. At the same time, it reveals that the steep decay and 
the following X-ray re-brightening have completely differ- 
ent spectroscopic properties (Fig. [5}: this, together with the 
temporal behaviour, strongly suggests that we actually see 
two different emission components ov erlapping for a sma ll 
time interval, as was first suggested bv lNousek et al.l (2006). 

The small overlap in time of the two components is the 
key ingredient that observationally allows the detection of 
the rising phase: this can be produced by either a steeper 



4 S ee the Swift-XRT l ight-curve repository, lEvans et ah! (120091 ') 
and lEvans et al.l 11200711 . 



to 


a 


X 2 /dof 


Model 


0(b) 


-7.1 ± 0.7 


1.7/4 




109 ± 89 (s) 


-4.2 ± 2.4 


2.5/5 




154 ± 13 (s) 


-3 


3.2/4 


Adiabatic cooling 


200 (s) 


-1 


42.1/4 


High latitude emission 



Table 6. Best fit parameters and statistical information for a 
simple power-law fit to the E p decay with time starting from 200 
s after trigger: E p ex. (t — to) a . 



5 DISCUSSION 

In GRB 081028 we have the unique opportunity to observe a 
smoothly rising X-ray afterglow after the steep decay: this is 
the first (and unique up to July 2009) long GRB Swift-XRT 
light-curve where a ris e with completely different properties 
to ty pical X-ray flares |Chincarini et al1l2007l . iFalcone et al.l 
|2007| ) is seen at t > 10 ks. At this e poch, canonical X-ray 
light-curves (e.g.. iNousek et al"1l2006l ) typically show a shal- 
low decay behaviour with flares superimposed in a few cases 
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Figure 10. Qualitative description of the spectral evolution with 
time detected in GRB 081028 from the prompt to the steep decay 
phase: the peak energy (-Bp) moves to lower energies while both 
the high and low energy components soften with time. Arbitrary 
flux density units are used. 



than usual steep decay or a delayed onset of the second com- 
ponent. We tested both possibilities comparing GRB 081028 
properties against a sample of 32 XRT light-curves of GRBs 
with known redshift and for which the steep-flat-normal de- 
cay transitions can be easily identified. While 63% of the 
GRBs are steeper than GRB 081028 (on ~ 2), no GRB 
in the sample shows a rest frame steep-to-flat transition 
time greater than 1 ks, confirming in this way the "delayed- 
second-component" scenario. Alternatively, the peculiarity 
of GRB 081028 could reside in a steeper than usual rise of 
the second component: unfortunately this possibility cannot 
be tested. 

This section is organised as follows: in Sect. 15. ll we dis- 
cuss the spectral evolution during the prompt and steep de- 
cay phases in the context of different interpretations. The 
afterglow modelling of Sect. 15.21 favours an off-axis geome- 
try: however, this seems to suggest a different physical origin 
of the prompt plus steep decay and late re-brightening com- 
ponents. This topic is further investigated from the prompt 
efficiency perspective in Sect. 15.31 



5.1 Spectral evolution during the prompt and 
steep decay emission 

The evolution of the peak energy E p of the vF v spectrum 
from the 7-ray to the X-ray band described in Sect. l4.6l offers 
the opportunity to constrain the mechanism responsible for 
the steep decay emission. 

Spectral evolution through the prompt and steep de- 
cay phase has been noted previously, with the E p track- 
ing both the overall burst behaviour and in dividual prompt 
pulse structures (see e.g.. iPeng et al. I l2009l for a recent time 
resolved spectra l anal ysis of prompt pulses). In particular, 
Yonetoku et~afl i|2008|) find E p oc t~~ 3 for GRB060904A; 
Mangano et alj (|2007l ) model the prompt to steep decay 
transition of GRB 060614 with a Band (or cut-off power-law) 
spectr al mo del with E D evolvin g as t~~ 2 , while iGodet et all 
l|2007t l and iGoad et all l|2007bl ) report on the evolution of 
the E p through the XRT energy band during single X-ray 
flares in GRB 050822 and GRB 051117, respectively. A de- 



caying E p was a lso observed during the 0.3-10 keV emission 
of GRB 070616 (|Starling et alj|2008l) . 

The detection of strong spectral evolution violates the 
prediction of the curvature effect in its s imple st formula- 
tion as found by IZhang. Liang fc Zhand ((2007) in 75% of 
the analysed GRBs tails: this model assumes the instanta- 
neous spectrum at the end of the prompt emission to be 
a simple power-law of spectral index j3 and predicts the 
a = 2 + (5 relation, where j3 is not supposed to var y (see 
e.g.. lFenimore et al.|[l996l : iKumar fc Panaitescull2000r i. The 
curvature effect of a comoving Band spectrum predicts in- 
stead Ep oc t -1 and a time dependent a = 2 + P rela- 
tion (see e.g.. lGenet fc Granotll2009l ; IZhang, Liang fc Wand 
120091 ): from Fig. [3 lower panel and Table El it is apparent 
that the observed E p oc t~ 7 ' 1±0,7 is inconsistent with the 
predicted behaviour even when we force the zero time of the 
power- law fit model to be to = 200 s, peak time of the last 
pu lse detected in the 15-150 keV energy range, as prescribed 
bv lLiang et alj l|2006l ). However, a more realistic version of 
the HLE might still fit the data: a detailed modelling is be- 
yond the scope of this paper and will be explored in a future 
work. 

The adiabatic expansion cooling of the gamma-ray pro- 
ducing source, which lies within an angle of I/7 (where 7 is 
the Lorentz factor of the fireball) to the observer line of sight, 
has also been recently proposed as a possible mechanism 
respo nsible for the steep decay l|Barniol Duran fc Kumarl 
2009). This process gives a faster temporal evolution of the 
break frequency as it passes through the X-ray band: typ- 
ically Ep oc t -3 . Two fits to the data have been done, the 
first fixing the break evolution to t~ 3 and the other one leav- 
ing to and the break temporal evolution as free parameters. 
Both fits are consistent with the adiabatic cooling expecta- 
tion and set to close to the beginning of the last pulse in 
the BAT light-curve (see Table EI- However, the adiabatic 
expansion cooling of a thin ejecta predicts a light-curve de- 
cay that is linked to the spectral index (5 by the relation 
a = 3/3 + 3, where a is the index of the power-law decay. 
Since a b s ~ 3.6 this would imply ~ 0.2 which is much 
harder than observed (Sect. |4~4|) . This makes the adiabatic 
cooling explanation unlikely. 

Both the curvature effect and the adiabatic model as- 
sume an abrupt switch-off of the source after the end of 
the prompt emission: the inconsistency of observations with 
both models argues against this conclusion and favours mod- 
els where the X-ray steep decay emission receives an im- 
portant contribution from the continuation of the central 
engine activity. In this case, the steep decay radiation re- 
flects (at least partially) the decrease in power of the GRB 
jet. An interesting possibility is given by a decrease of 
power originating from a decrease i n the mass accretion rate 
l|Kumar. Naravan fc Johnsonll2008l ). 

Alternatively, the observed spectral softening could be 
caused by cooling of the plasma whose cooling frequency 
identified with E D dec reases with time as suggested by 
IZhang. Liang fc Zhand l|2007h . 

While the spectral peak is moving, we also observe a 
softening of the spectrum at frequencies both below and 
above the peak when our data allow us to constrain the low 
and high energy slopes of a comoving Band spectrum. A 
softening of the low energy index in addition to the E p evolu- 
tion has been already observed in the combined BAT+XRT 
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analysis of GRB 070616 (jStarling et al.ll2008l . their Fig. 5). 
This result is consistent with the finding that while short 
GRBs have a low energy spectral component harder than 
long GRBs (i.e., |a B , short | < |Q B ,iong|, where ae is the low 
energy photon index of the iBandetal .11 19931 function), no 
difference is found in the ob distribution of the two classes 
of GRBs when onl y the first 1-2 s of long GRB prompt emis- 
sion is considered ( Ghirlanda et al.1l2009l ): a soft evolution of 
the ob parameter with time during the 7-ray prompt emis- 
sion of long GRBs is therefore required. Our analysis extends 
this result to the X-ray regime and indicates the softening of 
both the high and low spectral components from the prompt 
to the steep decay phase. The overall spectral evolution is 
qualitatively represented in Fig. 1101 



5.2 Afterglow modelling 

5.2.1 Failure of the dust scattering, reverse shock and 
onset of the afterglow models 

This subsection is devoted to the analysis of the X-ray re- 
brightening in the framework of a number of different theo- 
ries put forward to explain the shallow decay phase of GRB 
afterglows. 

According to the dust scattering model (|Shao fc Pail 

l2007h the shallow phase is due to prompt photons scat- 
tered by dust grains in the burst surroundings: this models 
predicts a strong spectral softening with time and a non- 
negligible amount of dust extinction which are usually not 
observed (Shcn et al ] |2009h . Both predictions are inconsis- 
tent with our data. 

A spherical flow is expected to give rise to a peak of 
emission when the s pectral peak enters t he energy band of 
observation (see e.g.. lGranot fc SarHl2002r ): the SED analysis 
of Sec. 14.51 clearly shows that E p was already below the X- 
ray band during the X-ray rising phase, well before the peak, 
thus ruling out the passage of the break frequency through 
the X-ray band as an explanation of the peak in the X-ray 
ligh t-curve. 

ISari fc Piranl (|l999T ) argue that the reverse shock has a 
much lower temperature and is consequently expected to ra- 
diate at lower frequencies than the forward shock, even if it 
contains an amount of energy comparable to the GRB itself, 
making a reverse shock orig in of the X-ray re-brightening 
unlike ly. However, following iGenet. Daigne fc Mochkovitchl 
|2007l ). in the case of ejecta having a tail of Lorentz factor 
decreasing to low values, if a large amount of the energy dis- 
sipated in the shock (e e near its equipartition value) is trans- 
ferred to only a fraction of electrons (typically £ e ~ 10 -2 ), 
then the reverse shock radiates in X-rays. In this case, it can 
also produce a plateau or re-brightening, the latter being 
more often obtained in a constant density external medium, 
that qualitatively agrees with the GRB 081028 afterglow. 

Alternatively, the detected light-curve peak could be 
the onset of the afterglow: in this scenario, the rising (de- 
caying) flux is to be interpreted as pre-deceleration (post- 
deceleration) forward shock synchrotron emission. The ob- 
served break frequency scaling oc t~ 2 - 6±0 - 2 ; s inconsistent 
with the expected cooling frequency evolution u c oc t~ x ^ 2 or 
v c oc t 1 / 2 for an ISM or a w ind environment, respectively 
(see e.g. iGranot fc Sarii r2002). We therefore consider a fast 
cooling scenario where Vb = ^m- The initial afterglow sig- 



nal from a thick shell is likely to ov erlap in time with the 
prompt emission (|Sari fc Piranlll999h , so that it would have 
been difficult to see the smoothly rising X-ray re-brightening 
of GRB 081028. For this reason only the onset of the for- 
ward s hock produced by t hin shells will be discussed. Fol- 
lowing ISari fc Piranl (|l999T l , the observed peak of the X-ray 
re-brightening implies a low initial fireball Lorentz factor 
70 ~ 75(noe 7 ,o.2) -1//8 , where no = n/(l cm -3 ) is the circum- 
burst medium density and e 7 ,o.2 = e 7 /0.2 is the radiative 
efficiency. Since the X-ray frequencies are always above the 
injection frequency v m , the X-ray light-curve should be pro- 
portional to t 2 j(t) 4+2p : during the pre-deceleration phase 
this means Fx oc t 2 for an ISM and Fx oc t° for a wind. 
The ISM scaling is consistent with the observed power-law 
scaling oc i 18±0 3 if a sharp transition between the rising 
and the decaying part of the re-brightening is required. The 
asymptotic value of the power-law index during the rising 
phase is instead steeper than 2, as indicated by the fit of the 
re-brightening where the smoothing parameter is left free 
to vary: oc £ 4,5±3 ' 3 (see Table [3] for details). The injection 
frequency is expected to scale as u m oc 7 (£) 4 - fc t- fc/2 , where 
the density profile scales as R~ k . This implies that for radii 
R < Rj (or t < tj) v m oc t° for an ISM and v m oc t^ 1 for 
a wind, while for R > R 1 (t > t 7 ) the fireball experiences a 
self-similar deceleration phase where 7 oc t~ 3 ^ 8 for an ISM 
and 7 oc i _1//4 for a wind, and i/ m oc t _3//2 in both cases. _R 7 
is the radius where a surrounding mass smaller than the shell 
rest frame mass by a factor 70 has been swept up; t 1 is the 
corresponding time: for GRB 081028 t 1 ~ 20 ks (observed 
peak of the re-brightening). While for t > t 7 the observed 
evolution of the break frequency is marginally consistent 
with i -3 / 2 , it is hard to reconcile the observed v m oc t~ a 
with a ~ 2.6 — 2.4 decay with the expected constant be- 
haviour or oc t -1 decay for t < t- 1 . This argument makes the 
interpretation of the re-brightening as onset of the forward 
shock somewhat contrived. Moreover, the identification of 
t = 20 ks with the deceleration time is also disfavoured by 
the earlier very flat optical light-curve. An alternative ex- 
planation is discussed in the next subsection. 



5.2.2 The off-axis scenario 

For a simple model of a point source at an angle of 6 from 
the line of sight, moving at a Lorentz factor 7 > 1 with 
7 oc R~ m l 2 , where R is its radius, the observed time is 
given by: 



t = 



1 



R 



2c7 2 \ 1 + m 



+ 1 



(4) 



The peak in the light curve occurs when the beaming cone 
widens enough to engulf the line of sight, 7(t pea k) ~ 1/0, 
so that before the peak t w R8 2 /2c oc R. We consider an 
external density that scales as R~ k (with k < 4) for which 
m = 3 — k. When the line of sight is outside the jet aper- 
ture, at an angle 6 from the outer edge of the jet, the emis- 
sion can be approximated to zeroth order as arising from 
a point source loca ted at an angle 9 from the line of sight 
l|Granot et al.ll2002l l. We have: 



to 
t 



v 



1 — j3 cos ( 



«aft 



1+7 2 



(•») 
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where (5 ■ 



-2U/ 2 



v/c and the subscript indicates 



the 8 = (on-axis) condition. The observed flux is given by 



F v {6,t) « at tt F v/a (0,at) 



(6) 



and peaks when 7 ~ 1/8. In the following we use the no- 
tations a a f t ~ 1/(1 + 7 2 # 2 ); a for the particular case where 
7 = To (where To is the initial Lorentz factor of the fireball): 
a^l/{l+Tl8 2 ). 

For t <C tpoak, 7$ 2> 1 and therefore a a f t ~ (7$) 2 oc 
7~ 2 oc R 3 ~ k oc f 3_fe . In this condition the local emission 
from a spherically expanding shell and a jet would be rather 
similar to each other, and the usual scalings can be used for 
an on-axis viewing angle (e.g.. [Granot fc Sarill2002T ): 



^ m ,o oc ir s(4 - fc)/2 oc r s ' 2 

v c , oc J?( 3fc " 4 )/ 2 oc t (a*-*)/C«-») 
with respective off-axis frequencies: 

~ D (fc~6)/2 ,(fc-6)/2 



a Vc oc R 



(2+fc)/2 



OC t 



(2+fe)/2 



(T) 
(8) 

(9) 
(10) 



For t > tpcak, a a ft ~ 1 and v « z^o, so that the break fre- 
quencies have their familiar temporal scaling for a spherical 
flow (eq. [7] and 00. 

For a uniform external medium (fc = 0), v c oc t and 
i _1//2 before and after the peak, respectively, while for a 
stellar wind environment (fc = 2) the corresponding tem- 
poral scalings are t 2 and t 1 ^ 2 . In both cases this is incon- 
sistent with the observed rapid decrease in the value of the 
break frequency (fb oc t~ 2,6 ) unless we require a very sharp 
increase in the magnetic field within the emitting region 
due to a large and sha rp increase in the external density 
l|Nakar fc Granod I2007T ). We consider this possibility un- 
likely (see Sect. 15.2. lj> . 

Alternatively, the break frequency could be v m , for a 
fast cooling spectrum where v c is both below v m and below 
the optical. In this case, for t < t pca k we have v m oc t~ 3 
(t~ 2 ) for afc = (fc = 2) environment; after the peak v„ 



t 



-3/2 



independent of fc. Since we observe Ub oc t 



(or 



Vb oc (t — to)~ 2 ' :i±lhl with to = 2ks) over about a decade in 
time around the light-curve peak, this is consistent with the 
expectations for a reasonable value of fc. 

Constraints on the model parameters are derived as 
follows: given that we see only one break frequency in 
our SEDs, which we identify with i/ m , we must require 
v c < v op t{~ 10 15 Hz). The tightest constraints are derived 
at tp Ca k, when u c reaches its maximum value (it increases 
with time before t PC ak and decrease s with time after i pC ak 
for fc < 4/3). From IGranot fc Saril \200t \. their Table 2, 
spectral break 11, this means: 



4 /2 "o^ 4 (l + F) 2 >10- 



pl/2 



(11) 



where eb is the fraction of the downstream (within the 
shocked region) internal energy going into the magnetic 
field; no = n/(l cm -3 ) is the external medium density; 

5 While these expressions are derived for a spherical flow, they are 
reasonably valid even after the jet break time tj et as long as there 
is relatively very littl e lateral expansion as shown by numerical 
simu lations (see e.g., Granot ct al. 2001c; Zhang & Mac FadvenI 
2009 and references therein). 



Ek,54, = -Ek,iso/(10 54 ergs) is the isotropic kinetic energy; 

Y is the Compton parameter which for fast coo ling reads 

Y w [(1 +4e e /e s ) 1/2 - l]/2, i|Sari fc Esinl l200lh : e c is the 
fraction of the internal energy that is given just behind the 
shock front to relativistic electrons that form a power-law 
distribution of energies: 7V C oc r y~ p for 7 ma x > 7c > 7mm- 
Assuming equipartition (e e = es = 1/3), Y ~ 0.62, Eg. 1111 
translates into: 



n > 2 x IQ~ S E 



'1/2 
fc,54 



(12) 



For an efficiency of conversion of the kinetic to gamma-rays 
energy e 7 = 1% the observed _E 7l i so = 1.1 x 10 53 erg (see 
Sect. g^J implies: n > 6 x 10" 4 . 

Using the best fit simple power-law models for the 
break frequency evolution with tim e of Sect. | 4.5I we hav e 
z^b(H2ks) ~ 1.5 x 10 14 Hz. Following IGranot fc Saril i|2002l ). 
their Table 2, spectral break 9, this means (a value that 
roughly agrees with the results for a range of values for p 
derived below is adopted): 



,1/2 



-1/2 



(13) 



where e e = £e7m/(7e) and £ e is the fraction of accelerated 
electrons. The value of p is p — 1.97±0.03 with intrinsic red- 
dening E(B - V) z = 0.03 (x 2 / dof = 135/138). Freezing the 
intrinsic reddening to E(B~V) Z = 0.06 gives p = 2.03±0.02 
(x 2 / dof = 140.8/139) while freezing it to E(B~V) Z = 0.08 
gives p = 2.08 ± 0.02 (x 2 / dof = 158.6/139). We thus take 
p = 2.0 ± 0.1. In particular, we calculate the range of values 
obtained for the microphysical parameters in the three cases 
p = 2.1, p = 2 and p = 1.9 since the expression of e e changes 
when p > 2, p = 2 and p < 2 l|Granot. Konigl fc PirarJ 
l2006t) : 



w(p-2)/(p-l) 

1/ ln(7 max /7min) 

(2-p)/(p-l)( 7n 



l/7n 



p > 2 
p=2 
p < 2 



(14) 



7max is obtained by equating the acceleration and cooling 
times of an electron, and y/3q e /(a T B'(l + Y)). 

Calculating the magnetic field value by B' = 
7aftC^/327reB nm p and assuming no = 1, e e = 0.3, 
as = 0.1 and 7a f t = 30 we obtain 

7max ^ 10 7 . Taking 
~ 500 (obtained for p ~ 2.1), (7 m i a / 7mM£ ) ~ 5 x 10 5 
(given the way this ratio appears in equation (I14|l - either 
in a logarithm or with a power 2 — p = 0.1 in our case 
- the dependence of the ratio e e /e e on it is very weak, 
and variations in its value have only a small effect). 
Then, since for p = 2.1, (p - 2)/(p - 1) ~ 0.1, and 
for p = 2, 1/ ln(7 ma x/ 7 min) ~ 0.1, for p > 2 we obtain 
(ee/Ce) 2 ^]/ 2 ~ 0.2. From the equipartition value - giving the 
maximum possible values e e /£e = es = 1/3 - we obtain an 
upper limit on the fraction of accelerated electrons: £ e < 0.3. 
For p = 1.9 we have (2 - p)/(p - l)(7min/7max) 2 ~ p ~ 0.04, 



and then (e e /£ e ) e 



21/2 



1.25, and then ^ < 0.2. The 



constraint on the microphysical parameters being very close 
in all cases, the exact value of p is then not of primary 
importance and the approximation p = 2.0 ±0.1 is then 
consistent. 

The evolution of the peak frequency being consis- 
tent with an off-axis interpretation of the afterglow, we 
further test this scenario by deriving the viewing and 
half-opening angle of the jet. The jet break time is 
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given by ISari Piran fc Halpernl <|l999h for the ISM and 
IChevalier fc Lil (|2000h for the wind environments: 



1.2 (1 + z) 



B.-4 



1/3 



0.1 



8/3 



days (fc = 0) 



6.25 (1 + z) (ff) (ff) 4 days (fc = 2) 



(15) 



From Table [3] we read a post-break power-law decay index 
6 = 2.1 ±0.1 (e= 2.3 ±0.1) if t ict ~ t pca k (tjet = *br„). Both 
are consistent with being post-jet break decay indices. We 
therefore conservatively assume tjet < 1 day, which leads to: 

-1/8 



A0 < < 



0.055 h^M. ra d 



(fc = 0) 



0.045 



-1/4 



(fc = 2) 



fif) rad 

Evaluating Eq. 9 of iNousek et all (|2006T ) at i 
7 ~ 1/9 we obtain: 

-1/8 



''peak 



7(tpeak) 



0.03 



0.03 



fir, i 
no 



rad 



-1/4 



rad 



(fc = 0) 



(ft = 2) 



(16) 



when 



(17) 



Using Eq. [12] for the ISM environment we finally have 9 > 
0.014£~ 3 { 16 rad. From the comparison of Eq.[l7]and Eq.fTol 
it is apparent that 9 > A9/2. Moreover, the slope of the 
rising part of the re-brightening of the afterglow is ~ 1.8, 
which is in rough agreement with the rising slope of t he re - 
brightening obtained from model 3 of iGranot et alj (|2002T ) 
- see their Fig. 2 - for 9 ~ 3A9. This is consistent with 
9 > A9/2. 

The off-axis interpretation implies that the value of the 
observed gamma-ray isotropic energy i? 7j j SOi fl corresponds 
to an actual on-axis input of E. 



-3i 



"y,iso,0 



'E. 



>7 ,iso,f3 



if e < 



10 53 erg, this may lead to very high energy output for this 
burst, which may be unphysical. It is therefore important to 
obtain limits on the Lorentz factor of the prompt emission, 
since a -1 ~ 1 + T 2 9 2 . Lower lim its to Fo can be obtained 
following ITithwick fc Sari |200ll) . requiring the medium to 
be optically thin to annihilation of photon pairs (Eq. I18[) 
and to scattering of photons by pair-created electrons and 
positrons (Eq. [T^lFI : 



-1/(2/33+2) 



J- mil 



1 150 keV \ 
\ m e c 2 J 



(fe-l)/(20B+2) 



n * (1 + z )(l-/3 B )/(/3 B + l) 

1/2 9 <A9 



/{ ( 0ji ) 1 /(23b+2) a -(/3 B +2)/(2(^ B +l)) #>A6> ^ 



<L/CS B +3) q + ^(/3b-1)/(/3b+3) 

-a/</»B-H» £) < A 



(19) 



where /3b is the high energy photon index of the prompt 

Band spe ctrum. 

From lfBlandford fc McKee] (|l976h . the Lorentz factor at 
the deceleration radius and at the peak of the re-brightening 

6 See Appendix [A] for a complete derivation of Eq. 1181 and 1191 



can be related by 7(i? pcak ) = 7(-Rdec)(i? P cak/-Rdec) (3 k)/2 - 
The Lorentz factor at the deceleration radius is a factor 
g < 1 of the Lorentz factor of the prompt emission To. 
Combining this with a -1 = 1 + Tq9 2 and 9 — l/7(i pea k), 
we obtain the following expression for the parameter a: 



R 



peak 



(20) 



Since g < 1/2, and Rdec < Rpcak, we have a -1 > 5 which, 
when substituted in equation [18] and [19] and keeping the 
strongest constraint, implies To > 46. To consider the other 
extreme case, where the deceleration time is ~ Tgrb, one 
should be careful in translating the ratio of radii to ratio 
of times: for a prompt emission with a single pulse, the 
duration of the GRB Tgrb is the duration of the pulse, 
which changes with the parameter a from on-axis to off- 
axis, as then does idee- We can therefore use off-axis values 
of the time t ~ R which means Rpe^/Rdec ~ t pca k/£dcc ~ 
£ P e a k/7GRB in our case here. Since i pca k ~ 2 x 10 4 s and 
Tgrb = 264.3 s (we identify the duration of the GRB with 
the T 90 parameter), a' 1 > 300 for fc = 2 (then T > 230) 
and a' 1 > 1.7 x 10 3 (To > 17 x 10 3 ) for k = 0. In the 
case of a prompt emission with several pulses, as it is the 
case for GRB081028, each pulse duration increases by a 
factor a -1 from on-axis to off-axis, however the total du- 
ration of the burst does not increase much, approximately 
by a factor of order unity, since the enlargement of pulses 
is somewhat cancelled by their overlapping. In this case, 
the GRB duration to consider is the on-axis one, for which 
t oc -R/7 2 oc R A ~ k \ since £ pca k is the limit between the on- 
axis and off-axis cases we can use t pca k oc Rp~ k k and then 

„ / t \ (3-fc)/(4-fc) 
^ = 1 + g- 2 (sgst) > 100 (or T > 136) for 

fc = and a" 1 > 36 (or T > 94) for fc = 2. 

The lower limit on the value of a -1 thus ranges be- 
tweerQ~ 5 and ~ 10 2 : this implies values of the isotropic on- 
axis gamma-ray energy output to range between -E 7 ,i so ,o ~ 
3 x 10 54 erg and £ 7 , iso ,o ~ 10 57 erg if 9 < A9 and even 
greater values for 9 > A9: between -E 7 ,i so ,o ~ 1.4 X 10 55 erg 
and -E 7 ,iso,o ~ 10 59 erg. These very high values could suggest 
that the observed prompt emission is from a different compo- 
nent than the observed afterglow emission. This possibility 
independently arises from the prompt efficiency study: the 
next section is dedicated to an investigation of this topic. 



5.3 Prompt efficiency 

The study of the efficiency of the conversion of the total 
initial energy into gamma-rays can in principle shed light 
on the physical mechanism at work. In the particular case of 
GRB 081028, this study helps us to understand if the prompt 
and afterglow emission originated from physically different 
regions. The first part of this sub-section is dedicated to the 
on-axis case; the second part to the off-axis case. 

Assuming that all energy not radiated in gamma-rays 
ends up in the kinetic energy E k of the afterglow, the im- 
portant parameters are the energy radiated in gamma-rays, 



7 Since GRB 081028 is composed of at least two pulses, we con- 
sider the most relevant case, when the observed off-axis duration 
of the prompt emission is close to the on-axis one. 
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_E 7 , the kinetic energy of t he afterglow, E k and a param- 
eter / = E k (10 hr)/g k , (|Granot. Konigl fc Piran] 120061 ; 
iFan fe Piranl 120061 . hereafter FP06) that accounts for en- 
ergy injection during the shallow decay phase (since energy 
injection is the most common explanation for this phase), 
where i?k,o is the initial kinetic energy of the afterglow, be- 
fore energy injection. Accounting for energy injection, the 
efficiency of the prompt emission reads: 



£ k ,o + E 1 1 + (/ - IK 



(21) 



where <f 7 ee E-,/(fE k , + E-,) = E-,/(E k (10 hr) + £ 7 ) is 
the prompt efficiency in the case of no energy injection. 
All the listed quantities are isotropic equivalent quanti- 
ties. The value of this parameter can be calculated with 
a good estimate of -Ek(10 hr) that can be obtained from 
the X-ray luminosity at 10 hours i f the X-ray frequency vx 
is ab ove both v m and v c (FP06; iLlovd-Ronning fc Zhan el 
|2004 hereafter LZ04). This is the case for GRB 081028 (see 
Sect. I5.2.2p which shows an isotropic X-ray luminosity of 
Lx,iso(10hr,obs) ~ (6.3 ± 1.0) x 10 47 ergs -1 . The calcula- 
tion of the kinetic energy is done following the prescriptions 
of FP06: unlike LZ04, they integrate their model over the 
observed energy band 0.3 — 10 keV and consider the effect 
of inverse Compton cooling. Equation (9) of FP06 gives the 
kinetic energy at ten hours: 



£ k (10hr)=J?4^ 6 +2 )(i±£ 



(2-p)/(p+2) 



(p-2)/(p+2) e 4(l-j>)/(p+2)^ 1 + y^4/(p+2) 



c B,-2 



(22) 



where .R = 9.2 x 10 52 [i(10hr)/T 90 ]T£r erg. This implies we 
need to make some assumptions on the microphysical pa- 
rameters e e , €b and Y. For the latter, as the afterglow is 
likely to be in fast cooling (see Sect. 15.272 i, then Y > 1 and 
we take Y ~ (e e /e s ) 1/2 following FP06. llVledvedevI (|2006h 
showed that during the prompt emission it is most likely 
that e e ~ ^/eB■ The values of the microphysical parame- 
ters being poorly constrained ("Sect. [5T2.2|) . we set e e = 0.3 
and e_e = 0.1, which is consistent with the values obtained 
in subsection 15.2.21 (see eq. Q~3] when £ e < 1 and eq. [14] 
and the paragraph below it). Taking p ~ 2, we thus ob- 
tain E k (10 hr) = 1.3 x 10 55 erg. Combined with the ob- 
served isotropic gamma-ray energy of the prompt emission 
E-y — 1.1 x 10 53 ergs, we have e 7 = 8.6 x 10 -3 (corresponding 
to a ratio £"k(10 hr)/E-, « 116): this is low, even compared 
to the values obtained by FP06 (their values being between 
0.89 and 0.01 - see their Table 1), which are already lower 
than previous estimates by LZ04. Now returning to the effi- 
ciency including energy injection, we can obtain an estimate 
of .Ek.o by using the previous formula but at the peak of the 
re-brightening and taking R — 9.2 x 10 52 erg (thus ignoring 
energy radiative losses since the end of the prompt emis- 
sion), which with its peak luminosity L pcak = 1.2 x 10 48 erg 
s _1 gives an initial kinetic energy injected into the afterglow 
£k,o = 1-16 x 10 55 erg and then an efficiency of the prompt 



emission which is as low as 



9.4 x 10" 3 . This calcula- 



tion assumes an on-axis geometry and accounts for energy 
injection. 

To strengthen the result of the above paragraph that 
the efficiency of GRB 081028 when considered on-axis is low, 
we analyse its prompt and afterglow fluencies and compare 



m 
or 

o 



E 2 




1 2 3 4 5 

Log[E 7 /t*L x . so (t)] 

Figure 11. Distribution of E 1 /tLx(t) with t = 10 hr rest frame, 
for the sample of 31 long GRBs detected by Swift with E 7 pro- 
vided by Amati et al. (2008). Black solid line: Gaussian best fit 
to the distribution. The dashed black line marks the position of 
GRB 081028 in the distribution, while the black solid arrow is 
pointed to the direction of increase of the radiative efficiency pa- 
rameter e 7 . 



them to a sample of Swift bursts from lZhang et aD (|2007bh . 
since fluences require no assumptions to be obtained. The 
prompt 1 - 10 4 keV gamma-ray fluencejfl of GRB 081028 is 
S-y ~ 8 x 10 _6 ergcm~ 2 and its afterglow X-ray fluence, 
calculated by 5*x ~ t peaik F v (tp eak ) to be consistent with 
IZhang et all (|2007bl) method, is Sx « 3 x 10" 7 ergcm" 2 , so 
that their ratio is S~,/S x w 26.7, plac i ng GR B 081028 in the 
lower part of Fig. 6 of IZhang et al.1 (|2007bl) . Compared to 
their sample of 31 Swift bursts, the 15 — 150 keV fluence of 
GRB081028, which is 3.2 x 10" 6 ergcm" 2 is well within their 
range of values (spanning from Sx.min ~ 8 x 10 _8 ergcm -2 
to Sx,max ~ 1.5 x 10 _5 ergcm -2 ; sixth column of their ta- 
ble 1), whereas its X-ray fluence is higher than most of them 
(see columns 6-9 of their table 2). It thus means that whereas 
GRB 081028 released as much energy in its prompt emission 
as most bursts, more kinetic energy was injected in its out- 
flow. This g ives a lower e fficienc y than most of the GRBs 
analysed by IZhang et all (|2007bh . consistent with the sce- 
nario above. Figure QT] clearly shows that this is likely to be 
extended to other Swtft long GRBs: at late afterglow epoch 
the X-ray band is above the cooling frequency and the X-ray 
luminosity is a good probe of the kinetic energy. In partic- 
ular E k oc Lx.iso (see Eq. 1221) : this means that high (low) 
values of the ratio i? 7 /L Xi i so are linked to high (low) values 
of radiative efficiency. 

The afterglow modelling of the previous section favours 
an off-axis geometry. In this case, considering that -E k ,i S o ~ 
10 55 erg, for the lower limit a" 1 ~ 5 (see Sect. I5.2.2p the 
efficiency of the prompt emission becomes e 7 ~ 0.23, which 
is a more usual value (it is in the middle of the efficiency dis- 
tribution of FP06). However, the upper limit of the range 



Depending on the high energy slope of the Band spectrum, 



we have S 7 ~ 6.6 X 10 erg cm 
9.5 X 10~ 6 erg cm" 2 for /3 B = -2.1. 



for fa 



-2.5 and S-, 
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of values for a gives an efficiency of 99% (when 9 < A9, 
and thus an even higher value for 9 > A9), which is ex- 
ceptionally high and very hard to reconcile with models of 
the prompt emission. This would suggest that the observed 
prompt emission is from a different component than the ob- 
served afterglow emission. 

An alternative way of achieving a more reasonable 
gamma-ray efficiency is if the observed prompt gamma-ray 
emission is from material along our line of sight, which has 
Sk.iso ~ S 7i i so , while the peak in the X-ray and optical light- 
curves at ~ 2 x 10 4 s is from a narrow jet-component pointed 
away from us that has a significantly higher £?k,iso- In this 
picture the afterglow emission of this material along our line 
of sight (and possibly also between our line of sight and the 
core of the off-axis jet component) could account for the very 
flat (almost constant flux) early optical emission (from the 
white light detection at 275 s, through the _R-band detection 
at f780 s, and the I-band detections at several thousand sec- 
onds). This early optical emission appears to be from a dif- 
ferent origin than the contemporaneous X-ray emission, and 
is most likely afterglow emission, regardless of the origin of 
the prompt emission: the observed X-ray and optical emis- 
sion in the time interval 1.8 ks ^ t ^ 9.5 ks implies a spectral 
index |/3ox| < 0.5. Conversely, assuming Pox = 0.5, the ex- 
pected X-ray contribution of the on-axis component at these 
times is w 3 x 10~ 4 mjy which is lower than the observed 
X-ray flux for t < 9 ks and comparable to the observed one 
at t ~ 9ks. 



6 SUMMARY AND CONCLUSIONS 

The 0.3-10 keV X-ray emission of GRB 081028 consists of 
a flat phase up to ~ 300 s (the XRT is likely to have cap- 
tured the prompt emission in the X-ray energy band) fol- 
lowed by a steep decay with flares superimposed extending 
to ~ 7000 s (component 1). The light-curve then shows a 
re-brightening which starts to rise at t ~ 8000 s and peaks 
around 20 ks (component 2). The different spectral and tem- 
poral properties strongly characterise the XRT signal as due 
to two distinct emission components. However, their further 
characterisation as emission coming from physically distinct 
regions is model dependent. 

The strong hard-to-soft evolution characterising the 
prompt and steep decay phase of GRB 081028 from trigger 
time to 1000 s is well modelled by a shifting Band function: 
the spectral peak energy evolves to lower values, decaying 
as £ poak oc t- 7 - 1±0 - 7 or £ poak oc (t - t Q )- i/2±2A when the 
zero-time of the power-law is allowed to vary: the best fit 
constrains this parameter to be to = 109 ± 89 s. In either 
case our results are not consistent with the oc t~ l behaviour 
predicted by the HLE in its simplest formulation. While a 
more realistic version of this model might still account for 
the observed £ pca k evolution, other possibilities must be in- 
vestigated as well: the adiabatic expansion cooling of the 
7-ray source predicts a steeper than observed light-curve 
decay and is therefore unlikely. While the peak is moving, 
a softening of both the low and high-energy portions of the 
spectrum is clearly detected. The failure of both the cur- 
vature effect and the adiabatic cooling argues against the 
abrupt switch-off of the GRB source after the prompt emis- 
sion and suggest the continuation of the central engine ac- 



tivity during the steep decay. An off-axis explanation may 
reconcile the high latitude emission or the adiabatic expan- 
sion cooling models with the data. This will be explored in 
a future work. 

GRB 081028 has afforded us the unprecedented oppor- 
tunity to track a smoothly rising X-ray afterglow after the 
steep decay: the rising phase of the emission component 
later accounting for the shallow light-curve phase is usually 
missed, being hidden by the steep decay which is the tail of 
the prompt emission both from the spectral and from the 
temporal point of view. The peculiarity of GRB 081028 lies 
in a small overlap in time between the steep decay and the 
following re-brightening caused by an unusual delay of the 
onset of the second component of emission. Contemporane- 
ous optical data allow the evolution of the SED during the 
re-brightening to be constrained: the spectral distribution is 
found to be best described by a photo-electrically absorbed 
smoothly broken power-law with a break frequency evolving 
from 1.6 x 10 15 Hz downward to the optical band. The break 
frequency can be identified with the injection frequency of 
a synchrotron spectrum in the fast cooling regime evolving 
as v\, oc t~ 2 - 6±0 - 2 . The intrinsic optical absorption is found 
to satisfy A v , z < 0.22. 

The observed break frequency scaling is inconsistent 
with the standard predictions of the onset of the forward 
shock emission even if this model is able to account for the 
temporal properties of the X-ray re-brightening (note that 
in this context the delay of the second emission component 
is due to a lower than usual fireball Lorentz factor or exter- 
nal medium density). Alternative scenarios have therefore 
been considered. While a dust scattering origin of the X-ray 
emission is ruled out since we lack observational evidence 
for a non-negligible dust extinction and strong spectral soft- 
ening, a reverse shock origin cannot be excluded. However, 
this can be accomplished only by requiring non-standard 
burst parameters: the ejecta should have a tail of Lorentz 
factors decreasing to low values; e e should be near equipar- 
tition; only a small fraction £ e ~ 10~ 2 of electrons should 
contribute to the emission. 

The predictions of the off-axis model have been dis- 
cussed in detail: according to this model a peak of emission 
is expected when the beaming cone widens enough to engulf 
the line of sight. The delayed onset of the second emission 
component is not a consequence of unusual intrinsic prop- 
erties of the GRB outflow but is instead an observational 
artifact, due to the off-axis condition. The observed evolu- 
tion of Vb is consistent with the expected evolution of the in- 
jection frequency of a fast cooling synchrotron spectrum for 
< k < 2. We interpret the light-curve properties as arising 
from an off-axis view, with 9 ~ 3A<9 and 9 ~ 0.03(-^)" 1/8 
for fc=0 (or 9 ~ 0.03(-§^)~ 1/4 for k=2), 9 being the angle 
from the outer edge of the jet and A9 the jet opening angle. 

In this scenario, the peculiarity of GRB 081028, or 
the reason why we do not observe more GRB 081028-likc 
events, may be attributed to the following reasons. Since 
GRB 081028 is a particularly bright (and therefore rare) 
event when viewed on-axis (with high on-axis .E iBO and L mo 
values), it is detectable by an off-axis observer even at the 
cosmological distance implied by its redshift z = 3.038. In 
addition, GRB 081028 appears to be characterized by a par- 
ticularly narrow jet, for which the ratio of the detectable off- 
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axis solid angle to on-axis solid angle is larger than for wider 
(but otherwise similar) jets. Finally, GRB 081028 might have 
a peculiar angular structure that is not representative of 
most GRBs, which would undermine the drawing of statis- 
tical conclusions under the assumption of a similar angular 
structure for most or all GRB jets. 

The radiative efficiency is one of the key parameters in 
GRB science: a precise estimate of this parameter would al- 
low one to distinguish between different models put forward 
to explain the observed emission. For the on-axis model, 
with e 7 ~ 10" 2 , the GRB 081028 efficiency turns out to be 
lower than the values obtained by FP06 and LZ04 for a sam- 
ple of pre- Swift GRBs: this directly implies that instead of 
having released as much energy in the prompt emission as 
most bursts of the two samples, GRB 081028 has a much 
greater kinetic energy injected in the outflow. Figure [11] 
clearly shows that this conclusion is likely to be extended 
to other Swift bursts with secure i? 7 ,iso measurement. This 
picture changes if we consider the off-axis interpretation: if 
the deceleration time is much longer than the prompt dura- 
tion the prompt and afterglow emission are consistent with 
originating from the same physical component and the effi- 
ciency of the burst is comparable to most bursts; if instead 
the deceleration time is close to the end of the prompt emis- 
sion, then the on-axis isotropic energy output would imply 
an extremely high efficiency of 99% which is very hard to ex- 
plain. This suggests that the prompt and afterglow emission 
come from different physical components. 

GRB 081028 demonstrates the evolution of GRB spec- 
tral properties from the onset of the explosion to ~ 10 6 s 
after trigger and shows that this is likely to be attributed to 
two distinctly contributing components of emission. These 
can be constrained only by prompt, broad-band coverage 
and good time resolution observations. 
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APPENDIX A: DERIVATION OF EQUATIONS 
18 AND 19 

As for the main article, the convention of a subscript (8) 
for on-axis (off-axis) quantities is used. A subscript " * " is 
added for the ca ses when 8 = Ad , Ad b eing the jet opening 
angle. Following iLithwick fc Sarll (|200ll ). their eq. 5 and 8, 
the lower limit to 70 due to photons annihilation reads: 



In 



l/(2/3 B +2) 



(1 +z)(1-<3b)/(/3b + 1) 



-En 



GS b -1)/(2/3b+2) 



(Al) 



while considering the scattering of photons by pair-created 
electrons and positrons: 



7» 



.l/(/3 B +3) 



(1 + 2 



,(/3 b -1)/(/3b+3) 



(A2) 



where: _E max = 150 keV for BAT observations; /3b is the 
high energy photon index of the pr ompt spectrum; z is th e 
redshift of the burst. From eq. 4 of ILithwick fc Saril (|200lh . 
the dimensionless quantity f can be re-written as: 

T0 = (2 - 1Xl ° } (JT /0.1s)(/fe-l) (A3) 

where: dt is the luminosity distance; 8T0 is the typical time 
scale of variability and /i,o is the on-axis number of pho- 
tons per second per square centimeter per MeV at the en- 
ergy of 1 MeV. The on-axis quantities must be now re- 
lated to the observed off-axis ones. In particular from eq. 
[5] directly follows STq = a 5Tg and vq = vg /a. The fluence 

^ = / ™f « EdN / dA while / = / « oc dN/dA. 
For a point source located at 8 > A8 dN/dA oc dQ oc 8 2 : 
this implies To = a~ 3 Tg, fa — a~ 2 fg. When 8 < A6 the size 
of the region significantly contributing to the observed emis- 
sion increases as 8 oc a" 1 : this translates into To = a~ 2 Tg, 
fo = a-7< 

continuity of the function at A8, it follows 



From the fact that t a ^ and requiring the 



to 



a Tg 



a* ' 



Tg 



8 < AO 
8 > A8 



where we remind the reader that a, = a(A8) 
TlA8 2 ). 

Substituting this result into eq. IA1I leads to: 



(A4) 



1/(1 + 



; l/(2/3 B +2) 



J- mm, 77 — 
,-1/2 



/ 150 kcV \ 
I me c 2 J 



(fe-l)/(2/3 B +2) 



(1 + z)< 1 -' 3 b)/(/3b+i) 



< A8 



(a,) 1/(2 ^ +2) a -CSB+2)/(2/3 B +2) Q >Ae 



(A5) 



Equation[l6]defines an upper limit to A8 that translates 
into a lower limit to a* considering that a(8) ~ /(l+j 2 8 2 ) ~ 
(76*) ~ 2 for -yd ^ 1. Inserting this information in the equation 
above and using STg = 70 s (variability time associated to 
the two pulses, Sect. l4~Tj) . d L = 17.4 Gpc, /3 B = 2.5, f& w 
6.8 x 10 6 , we finally obtain eq.QHand eq.[T9l 



APPENDIX B: TABLES 



This paper has been typeset from a TpX/ MpX file prepared 
by the author. 



c± = T, 



1/(/3b+3) 



a 
(a,) 



2/(/3 B 



(l + 2 ) CflB 

3) 



l/(/3 B +3) a -3/(/3 E 



-l)/(/3 B +3) 

6» < A8 

+3) 8 > A8 



(A6) 



The prompt spectrum of GRB 081028 does not allow to 
constrain the high energy photon index /3b, being consis- 
tent with a cut-off power-law (see Table [5]). Using /3b = 
— 2.5 (value we observe around 600 s, observer frame), 
/1 = 1.6 x 10~ 3 photons cm" 2 s _1 MeV" 1 . The observed 



evolution of /3b (see Sect. 14-6 P implies a harder high en- 
ergy spectrum at i < 600 s: using /3b = —2.1 we have 
fi = 3.6 x 10" 3 photons cm" 2 ^" 1 MeV" 1 . In the following 
/1 » 2 x 10" 3 photons cm" 2 s _1 MeV" 1 will be used. 
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Tmid 


Exp 


mag 


mag 






Flux 




Flux 












obs 


corr 






obs 




corr 








00 


(s) 










(mJy) 




(mJy) 








WHITE 


275.2 


147.4 


20.86 


20.70 


+0.46 


-0.32 


8.757 • 10~ 3 




1.023 ■ 10~ 2 


±3.030 


10" 


-3 


663.1 


19.4 


> 21.24 


> 20.87 






< 7.480 ■ 10" 


-3 


< 8.747 • 10" 3 








5174.2 


196.6 


> 21.02 


> 20.85 






< 7.619 ■ 10" 


-3 


< 8.910 ■ 10" 3 








6580.5 


139.6 


> 20.38 


> 20.21 






< 1.374 ■ 10" 


-2 


< 1.606 • 10~ 2 








101479.4 


8890.0 


22.94 


22.77 


+1.51 


-0.61 


1.301 • 10~ 3 




1.520 ■ 10~ 3 


±9.769 


10" 


-4 


124146.5 


8872.8 


21.46 


21.29 


+0.25 


-0.21 


5.078 • 10~ 3 




5.934 ■ 10" 3 


±1.056 


10" 


-3 


V 


185.9 


9.1 


> 16.99 


> 16.88 






< 5.811 ■ 10" 


-1 


< 6.430 • 10" 1 








366.6 


19.5 


> 18.60 


> 18.49 






< 1.319 ■ 10" 


-1 


< 1.460 • 10- 1 








712.7 


19.4 


> 18.06 


> 17.95 






< 2.169 ■ 10" 


-1 


< 2.400 • 10" 1 








4149.8 


196.6 


> 18.86 


> 18.75 






< 1.038 ■ 10" 


1 


< 1.149 • lO" 1 








5584.7 


196.6 


> 19.63 


> 19.52 






< 5.108 ■ 10" 


2 


< 5.6524 ■ 10~ 2 








11192.6 


598.5 


20.64 


20.54 


+0.53 


-0.35 


2.016 ■ 10~ 2 




2.225 ■ 10~ 2 


±7.752 


10" 


-3 


28542.1 


598.5 


19.82 


19.71 


+0.23 


-0.19 


4.304 • 10~ 2 




4.750 ■ 10~ 2 


±8.264 


10" 


-3 


45891.3 


598.6 


19.39 


19.28 


+0.16 


-0.14 


6.392 ■ 10~ 2 




7.054 ■ 10- 2 


±8.643 


10" 


-3 


57502.6 


598.6 


19.51 


19.40 


+0.22 


-0.18 


5.745 ■ 10- 2 




6.341 ■ 10~ 2 


±1.043 


10" 


-2 


101806.7 


8961.6 


> 19.94 


> 19.83 






< 3.839 ■ 10" 


-2 


< 4.249 • 10~ 2 








156450.1 


6146.6 


> 19.65 


> 19.55 






< 5.015 ■ 10' 


-2 


< 5.498 ■ 10~ 2 








B 


465.7 


19.4 


> 18.77 


> 18.63 






< 1.261 ■ 10" 


-1 


< 1.435 ■ 10" 1 








4969.6 


196.6 


> 21.54 


> 21.40 






< 9.837 ■ 10" 


-3 


< 1.119 • 10~ 2 








6404.7 


196.6 


> 20.66 


> 20.53 






< 2.212 ■ 10" 


-2 


< 2.493 ■ 10~ 2 








17796.8 


506.1 


20.34 


20.20 


+0.23 


-0.19 


2.971 • 10~ 2 




3.371 ■ 10~ 2 


±5.730 


10" 


-3 


35208.9 


483.4 


19.92 


19.78 


+0.18 


-0.15 


4.377- 10" 2 




4.966 ■ 10~ 2 


±6.686 


10" 


-3 


64072.7 


474.4 


20.80 


20.66 


+0.38 


-0.28 


1.943 ■ 10~ 2 




2.205 ■ 10~ 2 


±5.705 


10" 


-3 


101152.1 


8814.4 


> 21.54 


> 21.41 






< 9.837 ■ 10" 


-3 


< 1.10886 • 10~ 2 








155772.9 


6047.3 


> 22.70 


> 22.56 






< 3.380 ■ 10" 


-3 


< 3.845 ■ 10" 3 








U 


613.8 


19.5 


> 19.24 


> 19.07 






< 2.899 ■ 10" 


2 


< 3.390 • 10~ 2 








16976.9 


598.6 


> 20.58 


> 20.42 






< 8.438 ■ 10" 


-3 


< 9.778 • 10" 3 








23578.2 


511.4 


> 19.97 


> 19.80 






< 1.480 ■ 10" 


-2 


< 1.731 ■ 10- 2 








37243.3 


3516.9 


21.00 


20.83 


+0.38 


-0.28 


5.764 • 10~ 3 




6.727 ■ 10" 3 


±1.693 


10" 


-3 


66660.8 


3671.5 


> 20.63 


> 20.47 






< 8.058 ■ 10" 


-3 


< 9.337 • 10" 3 








123770.1 


8692.4 


> 21.08 


> 20.91 






< 5.324 ■ 10" 


-3 


< 6.226 ■ 10~ 3 








155516.4 


5699.7 


> 20.48 


> 20.31 






< 9.252 ■ 10" 


-3 


< 1.082 ■ 10~ 2 








UVW1 


416.4 


19.5 


> 22.03 


> 21.80 






< 1.412 ■ 10" 


-3 


< 1.746 • 10- 3 








589.5 


19.4 


> 18.68 


> 18.45 






< 3.090 ■ 10" 


-2 


< 3.819 • 10~ 2 








5994.8 


196.6 


> 22.63 


> 22.40 






< 8.128 ■ 10" 


-4 


< 1.005 • 10" 3 








33415.1 


885.6 


> 21.83 


> 21.60 






< 1.698 ■ 10" 


3 


< 2.099 ■ 10" 3 








40105.9 


885.6 


> 23.16 


> 22.93 






< 4.988 ■ 10" 


4 


< 6.165 • 10~ 4 








51671.4 


885.6 


> 22.86 


> 22.63 






< 6.576 ■ 10" 


4 


< 8.128 • 10~ 3 








65720.4 


4231.4 


> 22.19 


> 21.96 






< 1.219 ■ 10" 


3 


< 1.506 ■ 10" 3 








386961.8 


36503.7 


> 24.37 


> 24.14 






< 1.637 ■ 10" 


4 


< 2.023 • 10" 4 








733574.4 


42008.6 


> 23.92 


> 23.69 






< 2.477 ■ 10" 


4 


< 3.062 ■ 10~ 4 








UVM2 


564.5 


19.4 


> 19.44 


> 19.17 






< 1.477 ■ 10" 


-2 


< 2.663 • 10~ 3 








5789.7 


196.6 


> 21.30 


> 21.04 






< 2.663 ■ 10" 


-3 


< 3.383 ■ 10" 3 








29392.5 


771.3 


> 21.80 


> 21.53 






< 1.680 ■ 10" 


-3 


< 2.154 • 10~ 3 








54444.4 


4565.2 


> 23.48 


> 23.22 






< 3.576 ■ 10" 


4 


< 4.543 • 10" 4 








68159.5 


885.6 


> 21.30 


> 21.03 






< 2.663 ■ 10" 


3 


< 3.415 • 10" 3 
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Tmid 


Exp 


mag 


mag 


Flux 






Flux 










obs 


corr 


obs 






corr 






00 


(s) 






(mJy) 






(mJy) 






UVW2 


515.3 


19.5 


> 21.77 


> 21.47 


< 1.454 


10" 


-3 


< 1.917 


10" 


-3 


5380.0 


196.6 


> 21.48 


> 21.18 


< 1.900 


10" 


-3 


< 2.504 


10" 


-3 


56593.0 


885.6 


> 20.96 


> 20.65 


< 3.067 


10" 


-3 


< 4.080 


10" 


-3 _ 



Table Bl. Swift-WOT photometric set of GRB 081028. 3cr upper limits are provided in cases of non-detection. Column 1: observations 
mid-time since BAT trigger; column 2: exposure time; columns 3 and 5: observed magnitudes and fluxes; columns 4 and 8: extinction 
corrected magnitudes and fluxes; columns 6 and 7 report the errors on the extinction corrected magnitudes, while column 9 lists the 
errors on the extinction corrected flux. Only the Galactic extinction correction has been applied to the data. 



Tmid 


Filter 


Exp 


mag 


Flux 




mag 


Flux 












obs 


obs 




corr 


corr 






(s) 




(s) 




(mJy) 






(mJy) 






CrAO 


1779.84 


R 


23x60 


21.62 ±0.07 


(6.922 ± 0.446) ■ 10~ 


3 


21.545 ± 0.07 


(7.418 ± 0.478) 


io- 


3 


3585.60 


1 


30x60 


21.32 ± 0.09 


(7.560 ± 0.627) ■ 10~ 


3 


21.264 ± 0.09 


(7.961 ± 0.660) 


io- 


3 


5529.60 


1 


30x60 


21.43 ± 0.09 


(6.832 ± 0.566) ■ 10" 


3 


21.374 ± 0.09 


(7.193 ±0.596) 


io- 


3 


7473.60 


1 


30x60 


21.20 ±0.08 


(8.444 ± 0.622) ■ 10" 


3 


21.144 ± 0.08 


(8.444 ± 0.622) 


io- 


3 


9426.24 


I 


30x60 


20.66 ± 0.05 


(1.389 ± 0.064) ■ 10" 


2 


20.604 ± 0.05 


(1.462 ± 0.067) 


10- 


2 


GROND 


20880.0 


g' 




19.9 ±0.1 


(3.98 ±0.37) ■ 10~ 2 




19.79 ±0.1 


(4.406 ± 0.410) 


io- 


2 


20880.0 


r' 




19.3 ±0.1 


(6.92 ±0.64) -lO" 2 




19.22 ±0.1 


(7.454 ± 0.686) 


io- 


2 


20880.0 


i' 




19.2 ±0.1 


(7.59 ±0.70) ■ 10~ 2 




19.14 ±0.1 


(8.017 ±0.738) 


io- 


2 


20880.0 


t! 




19.1 ±0.1 


(8.38 ±0.77) ■ 10~ 2 




19.05 ±0.1 


(8.694 ±0.801) 


io- 


2 


20880.0 


J 




19.0 ±0.15 


(9.12 ± 1.26) ■ 10~ 2 




18.97 ±0.15 


(9.359 ± 1.293) 


io- 


2 


20880.0 


H 




18.7 ±0.15 


(1.202 ± 0.166) • 10" 


1 


18.68 ±0.15 


(1.221 ±0.169) 


io- 


1 


20880.0 


K 




19.0 ±0.15 


(9.12 ± 1.26) ■ 10~ 2 




19.00 ±0.15 


(9.135 ±0.502) 


io- 


2 


112680. 


g' 




21.26 ± 0.05 


(1.14 ±0.05) • 10~ 2 




21.15 ±0.05 


(1.259 ±0.058) 


io- 


2 


112680. 


r' 




20.49 ± 0.05 


(2.31 ±0.10) • 10~ 2 




20.41 ± 0.05 


(2.491 ±0.115) 


io- 


2 


112680. 


i' 




20.24 ± 0.05 


(2.91 ±0.13) • 10~ 2 




20.18 ±0.05 


(3.076 ±0.142) 


io- 


2 


112680. 


t! 




19.99 ± 0.05 


(3.66 ±0.17) • 10~ 2 




19.94 ±0.05 


(3.830 ±0.176) 


io- 


2 


112680. 


J 




19.6 ±0.1 


(5.25 ±0.48) • 10~ 2 




19.57 ±0.1 


(5.386 ± 0.496) 


10- 


2 


PAIR1TEL 


41133.2 


J 


1875.67 


17.78 ± 0.12 


(1.232 ± 0.126) ■ 10" 


1 


17.752 ± 0.12 


(1.264 ±0.140) 


io- 


1 


41133.2 


H 


1875.67 


16.91 ± 0.10 


(1.763 ± 0.162) ■ 10" 


1 


16.893 ± 0.10 


(1.791 ±0.165) 


io- 


1 


41133.2 


K s 


1875.67 


16.34 ±0.13 


(1.941 ± 0.232) ■ 10~ 


1 


16.3383 ±0.13 


(1.944 ±0.233) 


io- 


1 


44006.0 


J 


1844.28 


17.60 ±0.11 


(1.453 ± 0.147) ■ 10" 


1 


17.572 ± 0.11 


(1.492 ±0.151) 


io- 


1 


44006.0 


H 


1844.28 


16.83 ± 0.10 


(1.898 ± 0.174) ■ 10" 


1 


16.813 ± 0.10 


(1.928 ±0.178) 


io- 


1 


44006.0 


K s 


1844.28 


15.87 ±0.10 


(2.993 ± 0.276) ■ 10~ 


1 


15.8683 ±0.10 


(2.993 ± 0.276) 


10- 


1 


NOT 


19680. 


R 




19.23 ±0.03 


(6.255 ± 0.200) ■ 10~ 


2 


19.1545 ±0.03 


(6.706 ±0.185) 


10- 


2 



Table B2. Ground-based photometric set of GRB 081028. Column 1: observations mid-time since BAT trigger; column 2: photometric 
filter used; column 3: exposure; columns 4 and 5: observed magnitude and flux; columns 6 and 7: magnitudes and fluxes corrected for 
Galactic reddening. GROND data come from Clemens et al. (2008; 2008b). CrAO data come from Rumyantsev et al., (2008). 



